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1.0 SUMMARY 
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Significant differences exist in the noise generated by engines in 
flight and engines operating on the test stand. It has been observed 
that these differences are reduced by the use of an Inflow Control 
Structure in the static test configuration. It is the purpose of this 
contract (NASl-15085) to produce a design system for inflow control 
structures. This interim report describes the results of the develop- 
ment of the analytical models required for this design system. 

Characteristics of the fan inflow field are discussed based on observa- 
tions of fan blade mounted pressure transducer and inflow hot film 
data. Models were developed to study the effect of flow contraction and 
screening on inflow distortions, to describe the acoustic characteris- 
tics - transmission loss and directivity change - of honeycomb and per- 
forated plate, and to assess experimentally the directivity change due 
to a corner between honeycomb panels. 

The primary results of this study are: 

0 The distortions in the engine inflow field during static operation 
are due to atmospheric turbulence, the nacelle boundary layer, 
exhaust flow reingestion, and flow over stand structure, the 
ground plane and the engine casing. The dominant factors depend on 
the configuration and ambient conditions. 

0 For the JT9D engine on the PWA test stand, the dominant factor in 
the distortion field is flow over stand structure and ground 
effects. 

0 Flow contraction suppresses, initially, turbulent axial velocity 
distortions and magnifies turbulent transverse velocity distor- 
tions. 

0 Perforated plate and gauze screens preferentially suppress the 

axial components of both steady and turbulent velocity distortions 
to a degree determined by the screen pressure loss coefficient. 

0 Honeycomb screens preferentially suppress the transverse com- 
ponents of both steady and turbulent velocity distortions, to a 
degree determined by the length to diameter ratio of the honey- 
comb. Values of this ratio as small as 2 or 3 produce large reduc- 
tions in the transverse components. 

0 There is an observed tendency to isotropy in turbulent velocity 
fields. This relieves the anisotropy of the turbulent field found 
downstream of screens and limits the distortion of turbulence 
fields imposed by flow contraction. 



0 The acoustic transmission loss of perforated plate is controlled 
by the reactance of its acoustic impedance. A criterion is given 
for design of perforated plates having a transmission loss less 
than 1 dB in the frequency range of interest. 

0 The acoustic transmission loss of low solidity honeycomb screens 
is negligible for the external hemispherical design of inflow con- 
trol structures. 

0 A model for the direction change due to a corner between honeycomb 
panels is compared to measured data and the comparison is good. 

Based on these studies, certain qualitative design criteria have been 

identified. 


2 



2.0 INTRODUCTION 


Reduction of jet engine noise levels is a continuing evolutionary pro- 
cess which has resulted in significant noise reduction since the advent 
of the original commercial jet aircraft in the early 1950*s. Engine 
noise reduction features are usually verified by means of static or 
flight testing. Since the cost of flight testing is much greater than 
that of static testing, use of the latter technique results in signifi- 
cantly lower costs being accrued to the development of engine noise 
reduction features. As such, the use of static testing techniques 
should be fully exploited. The purpose of this contract is to develop, 
for turbofan engines, static testing techniques that result in measure- 
ments of fan noise that are representative of inflight levels. 

It has been noted by various observers that the noise produced by the 
fan of turbofan engines operating statically on the test stand is 
greater than that produced when the engine is operating under flight 
conditions. This may be a result of both fan blade passing tone and 
broadband levels being contaminated by extraneous noise sources present 
during static testing. As a consequence, predictions of flight noise 
levels using static data are consistently high. A no less important 
ramification is that flight noise sources cannot be identified from 
noise measurements from statically' operating engines. Correspondingly 
noise reduction techniques cannot be evaluated on the test stand in the 
presence of a contaminating noise floor. To identify the source of ex- 
traneous noise it is necessary to note that, statically, the fan inter- 
acts with a more distorted inflow field than it does when in flight. 
There are several distinct characteristics of the inflow field that are 
quite different when the engine is operating statically, each of which 
could produce extraneous noise. 

Firstly, the ambient atmospheric turbulence field in the vicinity of 
the ground is more turbulent than at higher altitudes. In addition, 
this turbulence field is convected through a very high flow contraction 
when the engine is operating statically, whereas, in flight, the turbu- 
lence field convects through a very small contraction on its way to the 
fan. This high flow contraction in the static case results in a distor- 
tion of the field in which the "eddies" are "stretched", then "chopped" 
by successive fan blades, producing "bursts" of discrete noise that are 
virtually absent in the flight operation of the engine. 

Secondly, in the static case, there exists stand structure and a ground 
plane that are sources of flow disturbances. Engine ingested air passes 
over the stand structure and the ground resulting in vortices and wakes 
in the inflow field. Engine case protuberances can also generate 
ingestible distortions. Usually these sources of distortion do not 
exist inflight since the air ingested into an engine has not passed 
over any boundaries. 
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Thirdly, the nacelle boundary layer is different during static opera- 
tions of the engine. In flight, at takeoff and approach conditions, air 
is drawn into the engine through a streamtube forward of the engine, 
which is slightly larger than the fan, causing the locus of the stagna- 
tion points to be near the leading edge of the nacelle. In static op- 
eration, however, air is ingested from a much wider streamtube causing 
the stagnation locus to move rearward on the outside of the flight na- 
celle. The acceleration of the flow from this stagnation locus around 
the lip region produces higher local velocities than that occuring in 
flight and consequently, a larger adverse pressure gradient along the 
nacelle inner wall than that which is encountered inflight. The bound- 
ary layer then becomes thicker in this region and may even separate 
(the use of bell mouth inlets during static testing is an attempt to 
alleviate this situation). Due to these differences in boundary layers, 
the area of the fan tip that is affected by the boundary layer during 
static operation is probably larger. The nacelle boundary layer in 
static operation is also affected by the presence of free stream flow 
distortions which can locally perturb it. It seems probable, then, that 
the nacelle boundary layer generated by static operation of the engine 
constitutes a larger flow disturbance than its flight counterpart and 
thus could be a cause of higher noise levels. 

Finally, since it is possible for the flow to be drawn from all angles 
by an engine on the test stand, turbulent flows from the jet plume may 
be reingested. This distortion source is not present in flight. These 
disturbances (i.e., atmospheric turbulence, ground plane and stand in- 
duced distortions, dissimilar nacelle boundary layer and jet plume re- 
ingestion) are considered to be the most important extraneous noise 
sources in static engine operation. In order to obtain useful static 
acoustic data, therefore, it is necessary to develop techniques which 
modify the inflow field so that, for all intents and purposes, the fan 
is operating as it would in flight. In the past, various techniques for 
accomplishing this simulation have been used, including mounting en- 
gines in wind tunnels and using devices upstream of the fan to condi- 
tion the inflow. Inflow Control Structures (ICS) for conditioning the 
flow have been mounted upstream of the engine by several investigators. 
A boundary layer suction system has also been used in conjunction with 
a screen. Both of these techniques have resulted in reduced radiated 
noise levels, indicating the reduction of inflow distortion. 

In view of the encouraging results achieved by the use of ICS's, the 
present contract was awarded for the purpose of developing an inflow 
control screen design procedure and a flight noise prediction procedure 
using data gathered from the static testing of engines equipped with 
such a structure. The use of wind tunnels, or other techniques for sim- 
ulating the flight environment are not considered under this contract. 
In addition, the problem of nacelle boundary layer modification is not 
addressed under this contract, but is certainly an important area for 
future work. The major components of the contract are: 
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Phase I Definition of Atmospheric Turbulence Characteristics and 
Engine Sensitivity Study. 

Phase II Development of Inflow Control Structure (ICS) Preliminary 
Design System. 

Phase III Interim Procedures Report Development and Coordination. 

This interim report describes the analytical models developed by the 
contractor that describe the physical processes involved in the passage 
of distortions and turbulence through contraction and screening and 
propagation of sound through screens. These models were developed under 
Phase II, Tasks A, C, D and E of the contract. 

The report contains a brief survey of the literature relevant to inflow 
control screens followed by a discussion of the various design consid- 
erations, both acoustic and aerodynamic. Next, in order that the upwash 
velocity field at the fan face may be quantified using blade mounted 
pressure transducer (BMT) data, the relationship between the unsteady 
velocity field at the fan face and the pressure field at a point on a 
blade is discussed. (The screen design system will be assessed using 
this information in a subsequent portion of the contract). The analyti- 
cal models of the effects of a screen and flow contraction on both 
steady and unsteady flow distortions are described. The acoustic char- 
acteristics of the screen components (including the transmission loss 
through the screen material and the directivity changes due to screen 
discontinuities) are then considered, followed by an experimental ver- 
ification of the directivity model. 
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3.0 SYMBOLS 


Pi Screen attenuation factor 
A. Defined in Section 10.0 
a Characteristic dimension 
Ouq Speed of sound 
B Screen process 
B Defined in Section 7.0 
b Wake width 
C Contraction process 

Contraction coefficient 
Drag coefficient 
Q Fresnel integral 
C Integration path 
Z) Length of contraction 
d Effective diameter of honeycomb 
d Distance between t-planes 
B Constant in Section 7.0 

Components of diagonal strain tensor 
F Velocity one dimensional power spectral density 
3* Fourier transform of velocity field 
■P Frequency 

G Geometric screen parameter 
H Hankel function 
7{ Defined in Section 8.0 
A Plate separation 
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3,0 Symbols (Cont'd) 

I Intensity 

Integrals 

J Defined in Appendix C 

K Screen pressure loss coefficient 

Defined in Section 8.0 
k Wavenumber 

L Length scale 

JL Thickness of honeycomb or screen 

JLf Contraction ratio 

Characteristic length in turbulence 
M Mach number 

rrt Defined in Appendix G 

N Filtering process 

n Index 

P Pressure power spectral density 


yD Acoustic pressure amplitude 

y> Instantaneous acoustic pressure 

Dynamic head 

Cyli ndrical coordinate system 
Pji, Reynolds number 

^ Resistive part of screen acoustic impedance 


f Position vector 

S Sears response function 

S Fresnel integral 


5 Screen sol i di ty 
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3.0 Symbols (Cont'd) 

T BMT response function 

T Defined in Section 10.0 
t Time 

H Mean velocity 

^ Acoustic particle velocity amplitude 

U Instaneous unsteady velocity 

)l Total velocity 

V Distortion velocity fiel d 

1/ Fan face distortion velocity field 

W Defined in Appendix F 

y Distance downstream of rod 

X Position vector 

^ Impedance 

z Complex number 

<x Screen flow angle ratio 

/3 Screen corner angle 

r" Three dimensional velocity spectrum tensor 

y Blade relative mean flow angle 

6 End correction 

c Small distance 

^ Phase difference 

Tj Lumley's RMS velocity ratio 
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3.0 Symbols (Cont'd) 

0 Angle 

'h! Post contraction wavenumber 

X Wavelength 

Mean squared velocity component ratios 
u Kinematic viscosity 

1 Distance 

p Air density 

a Open area ratio 

<f> Unsteady velocity potential 

X Reactive part of screen impedance 

SI' Stream function 

Q Total vorticity 

u Unsteady vorticity 

(p Radial frequency 


9 



Subscripts: 

^ Pre, contraction or screen 

<2. Atmospheric turbulence 

B Post, contraction or screen 

b Stand structure induced 

C ICS induced (5.0), Honeycomb cell (E) 

d Inlet induced 

F In flight 

I Incident 

/ Denotes velocity component 

^ Denotes wavenumber component 

/ Lower 

h Upwash 

R Reflected 

S Statically 

5 Screen 

~T Transmitted 

b Total 

u Upper 
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4.0 BACKGROUND AND LITERATURE REVIEW 

It has long been known that the flight and static test environments of 
an engine are different. For instance, in 1966, Sofrin and McCann (1) 
recognized the existence of noise generation mechanisms peculiar to 
ground operation of an engine. They identified "natural atmospheric air 
currents" and ground vortices as causes of inflow distortion that, upon 
interaction with a rotor, produce noise. They also surmised that the 
inflow into a fan in flight was not distortion free, but probably had a 
different distortion flow field than which existed during static 
testing. 

Subsequently, evidence accumulated (2, 3, 4, 5, 6) that showed not only 
were static and flight inflow distortions different, but that the sta- 
tic distortion field generated noise levels, especially blade passing 
tone levels (and possibly its harmonics) that were significantly 
greater than those generated in flight. However, while the source of 
the extraneous static noise is generally ascribed to the more disturbed 
inflow field, the particular distortion type responsible is not clearly 
defined. For example, Cumpsty and Lowrie (2) produced results that in- 
dicated that boundary layer changes were important at relative tip Mach 
numbers less than .85. In the flow outside the boundary layer, they 
suggest that distortions produced by the contraction of the ambient 
atmospheric turbulence field are important. Hanson's (7, 8, 9) measure- 
ments of inlet flow fields are explained by that author also as the 
result of contraction of the atmospheric turbulence, although, in some 
situations, he also identifies distortions due to flow over stand 
structure. Distortions have also been detected that were generated by 
protuberances on the exterior nacelle casing (10). Earlier work by 
Povinelli and Dittmar (11) identified distortions of the inflow by the 
installation structure. Hodder (12) examined among other things, the 
noise generated by a JT15D at two heights from the ground, and found 
the ground vortex interacting with the rotor to be a significant 
source. Recent P&WA/Boeing experience under the Joint Noise Reduction 
Program (JNRP) indicates that on a full size test stand, ground ef- 
fects, stand structure and possibly atmospheric turbulence contribute 
to inflow distortions. The above review indicates that any or all of 
the distortion sources can exist but to a varying degree depending on 
the details of the test stand configuration. 

In order to develop a design system for ICS's, it is necessary to quan- 
tify the effects of contraction and screening on inflow distortions 
that are convected into an engine inlet during static testing. Wind 
tunnel designers have used screens to suppress flow distortions, and as 
a result a body of literature has evolved describing the effects of 
screening on flow distortions. For example, Prandtl (16) had long ago 
determined that wind tunnel flows could be made steadier with gauze 
screens and honeycomb. He, and later Collar (17), and Batchelor (18) 
had attempted analyses of the phenomenon. However, it was Taylor and 
Batchelor (19) who derived a solution of small flow perturbations 
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convected through a gauze screen. They showed how two characteristics 
of the screen determine the degree of turbulence of steady distortion 
suppression. These are the pressure loss coefficient and the flow de- 
flection ratio (the amount of turning of the flow caused by the gauze). 
Honeycomb, as shown by Lumley (20) acts in a similar manner. A signifi- 
cant factor in the use of screens for flow distortion suppression in 
wind tunnels is the consideration of the turbulence and distortions 
generated by the suppressing screens themselves. To minimize self gen- 
erated turbulence and flow distortions, combinations of honeycomb and 
screening have been used (See for example, (21)), The body of litera- 
ture described in this paragraph was used as the basis for developing 
the models described herein, which account for screen effects on turbu- 
lence and steady distortions. 

The effects of contraction on turbulence and steady distortion have 
also been addressed in the literature, again motivated by the need for 
understanding these phenomena in the design of wind tunnels. Models for 
the effects of contraction were developed by Prandtl (16), Taylor (35), 
Kibner and Tucker (37), Batchelor and Proudman (38), and Pearson (39). 
These authors analyzed the problem using differing assumptions (See 
detailed discussion in section 7.2). Experimental results have also 
been reported in the literature, e.g., Townsend (42), Uberoi (38), 
Hussein and Ramjee (41) and Reynolds (43). Experimental results have 
been limited to contraction ratios that are small compared to those 
encountered in static engine operation. Comparisons of existing models 
with the available data show poor agreement. Postulated reasons for the 
discrepancies cited in the literature are discussed in section 7.2. In 
order to define contraction effects for the ICS design system required 
for this contract, the analytical results of the cited references will 
be used as a basis for the contraction model, and the results of the 
test program (Phase II, Task F) will be used to empirically correct the 
models. This interim report will describe the analytical development of 
the models. 

Most published work to date on the use of screens to better simulate 
fan inflow fields for fan tests has been performed using model scale 
fans. The ICS's tested to date have been performed using model scale 
fans. These ICS's have been of the in-duct type and the externally 
mounted hemispherical type. The in-duct screen has consisted of either 
a honeycomb panel mounted upstream of the rotor (12, 22, 23) or mul- 
tiple gauze screens at the inlet entrance (4) (It should be noted that 
the tests conducted by Jones et al was not a model test, but was con- 
ducted on a JT15D, a small turbofan, at a test stand). 

All except Leggat and Siddon observed a reduction in tone level, how- 
ever, these investigators did note a change in blade loading caused by 
the suppression of transverse velocities in a ground vortex. Suppres- 
sion of this vortex nearer its stagnation point resulted in noise re- 
ductions. In conjunction with the noise reduction, Hodder (12) observed 
a reduction in axial length scale, as determined by making auto 
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correlations of the hot wire signal obtained at the fan face. Hodder 
ascribed the noise reduction to a change in length scale caused by the 
screen. However, Jones, et al, (23) demonstrated that honeycomb screens 
mounted in the inlet duct have a substantial acoustic transmission loss 
associated with them, so there is some question as to the reason for 
the noise reductions when the screens are mounted in the inlet, i.e., 
are they due to source noise reduction, transmission loss across the 
screen, or a combination of both? 

The externally mounted hemispherical screens do not suffer from this 
problem of acoustic transmission loss to any great degree (27). In gen- 
eral, all screens of this type (13, 14, 23, 28, 29, 30) produced a sig- 
nificant reduction in tone noise, at subsonic relative tip Mach numbers, 
whether in the anechoic chamber, wind tunnel or on the test stand. The 
unsteadiness of the tone also decreased. These reductions were also ob- 
served with P&WA/Boeing JNRP screens mounted in front of a JT9D on the 
test stand. This noise reduction was accompanied by a reduction in tur- 
bulence level and scale. Only Shaw, et al , (13) show results in con- 
flict with the scale reduction, but in this case, the screen was tested 
in a rectangular wind tunnel that was not much larger in height than 
the screen diameter. This situation could result in distortions gener- 
ated by the screen itself. Based on this review of existing experience 
with inflow control structures, there appears to be a consensus that 
use of ICS's during static testing does result in a better simulation 
of flight fan noise generation mechanisms. 

While universal screen design criteria have not yet been determined, 
certain desirable characteristics have emerged based on the literature 
reviewed. These characteristics are summarized below: 

0 The minimization of acoustic transmission loss is enhanced by use 
of external ICS's as opposed to use of ICS's mounted in the engine 
inlet. 

0 Distortion suppression, primarily of transverse velocities, is 
effected by use of honeycomb. Honeycomb characteristics suggested 
by Ho, et al , (31) are that a length to diameter ratio of 5 to 10 
is necessary. However, Ginder (32) shows that increasing this 
ratio beyond 2 does not improve disturbance suppression. This 
latter view is born out, in some degree, by Prandtl (15). 

0 Axial velocity distortion suppression by use of gauzes or screens 
because of their relatively high pressure loss coefficients. 

0 Distortions generated by honeycomb structure can be suppressed by 
use of gauzes mounted downstream of the honeycomb, as suggested by 
the work in references ( 13 ), ( 21 ), ( 31 ) and ( 32 ). 

0 Gauzes, screens, and perforated plates can be mounted upstream of 
the honeycomb to protect it from damage by incoming debris. 
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Inflow control devices should be designed to avoid large struc- 
tural members and discontinuities (as noted by Lowrie and Newby 
(21) and PWA/Boeing experience) since they can give rise to self 
generated distortions entering the fan. 

As can be seen from the above review of the literature, the feasibility 
has been shown of using an inflow control structure during static test- 
ing to effect a better simulation of the fan inflight flow field. Cer- 
tain desirable characteristics of ICS designs have been identified, but 
as yet, a design procedure does not exist. Additionally, models des- 
cribing the effects of screening and contraction on steady wake and 
turbulence distortion exist, although verifying test data for contrac- 
tion ratios of the size encountered in static testing of fans does not 
exist. Additionally, there is only limited data on the effects of 
screening. Accordingly, the objective of the subject contract is to 
develop a procedure for design of ICS's by developing appropriate ana- 
lytical models of screening, contraction and acoustic transmission ef- 
fects and to provide the data required to empirically verify and/or 
correct these models. 
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5.0 INFLOW CONTROL STRUCTURE DESIGN CONSIDERATIONS 


The inflow control screen has to perform the function of suppression, 
in some degree, of ambient flow distortions. It has to fulfill this 
requirement in such a manner that it does not itself produce signifi- 
cant flow distortions and such that the acoustic data, gathered when 
using it, is meaningful. Furthermore, the inflow control screen should 
not significantly affect the engine operating line. The criteria speci- 
fying the suppression, generation and acoustic effects of the inflow 
control screen are examined in the following subsections. In the fol- 
lowing discussion, the term distortion velocity is used to indicate the 
fan face velocity induced beyond that which would exist if the fan were 
operating in isolation, in an undisturbed atmosphere. 

5.1 Distortion Suppression Criteria 

The use to which the acoustic data, gathered with an inflow control 
device in front of the fan is put, determines the criteria governing 
the design of the screen. Currently two uses are identifiable: 

a) The acoustic data are corrected to give corresponding flyover 
levels. 

b) The acoustic data are used to evaluate the effects of design 
changes on radiated noise levels. These changes may be made on 
full scale engine or on models. 

The first data use requires that all screens should produce a fan face 
velocity field that is substantially the same as that in flight. The 
second data use requires that the screen suppress the fan face distor- 
tion field to flight levels, or to as low a value as possible, depend- 
ing on whether the application is assessed of the impact of the change 
in flight noise levels or a scientific study where suppression of all 
inflow distortion sources is desired. To aid in considering what the 
criteria for simulation and minimization should be, a schematic diagram 
of the sequence of events involved in static and flight operation is 
shown in Figures 1 and 2. Considering Figure 1, four source types of 
fan inflow distortion are shown, i.e., atmospheric turbulence, distor- 
tions induced by stand structure, etc., distortions induced by the ICS 
and finally, inlet induced distortions. The boxes show the processes 
that alter these distortions, with C/aS for example indicating alter- 
ation of atmospheric turbulence by contraction from the atmosphere to 
the screen, in the static case. Bq. signifies its alteration by pas- 
sage through the screen and C^as signifies the alteration caused by 
further contraction from the screen to the fan face. The alterations of 
the other disturbance sources are indicated in similar fashion. The 
static case without an ICS is obtained by noting that there is no 
attenuation by the screen. The processess shown in Figure 1 for each 
source type can be written symbolically as shown below in Equations 1, 

2, 3 and 4. 
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Figure 2 


Sources of fan face distortion during flight operation^ 
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where the variables 14 s»l 4 s* '"eP'"esent distortion velocities as 

they exist at the source, represent the distortion ve- 
locities at the fan face under static conditions and the symbols 
C[] and S[] represent operators for alterations caused by 
contraction and screening. The total fan face distortion velocity is 
the summation of Equations 1 , 2, 3 and 4, i.e., 

~ ^/aS Kisj + 

^2cS ^cS 

Keference to Figure 2, which shows the sources of fan face distortion 
during flight operation, results in the following equation for the fan 
face distortion velocity 
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where Vp is the total fan face distortion velocity in flight, 

Vqf i VdF i^epresent distortion velocities as they exist at their 
source and, as before, the symbol C[J represents an operator for the 
alterations caused by contraction. 

Simplification of Equations 5 and 6 may be achieved by assuming that 
the self generated distortions VcF » and '^dF are zero. In 

flight this assumption is true away from the boundary layer, at small 
angles of attack, and for engines mounted not too close to the fuselage 
(Blade Mounted Transducer data from a JT9D engine mounted inboard on a 
747 verify these assumptions for this installation). Statically this 
assumption is true outside the boundary layer and then only if the in- 
let flow control screen does not generate distortions, a property that 
can be incorporated by proper design. This latter point is examined in 
the next subsection. 

To further refine the expressions in 5 and 6, if it is noted that the 
fluctuating fan blade forces are the dominant noise generation mechan- 
ism, then only that component of the external velocity field that con- 
tributes to the fluctuating blade forces, i.e. , the component of the 
distortion velocity perpendicular to the mean velocity, is of interest. 
In addition, this upwash velocity field can be considered to be made up 
of contributions from many distortions eddies of varying cross section 
and length. Only certain of these eddies contribute velocities that 
produce noise at blade passing frequency and its harmonics (although 
all unsteady components may produce broadband noise). The expressions 5 
and 6 can be rewritten to reflect these considerations, i.e.. 






where the subscript n denotes the upwash component of velocity and 
symbol //| } represents the operation of selecting only those 
"eddies" that give rise to blade passage frequency tone noise and its 
harmonics; i.e., it is mathematically equivalent to filtering in wave 
number space. 

Having evolvefJ Equations 7 and 8, the next step is to establish cri- 
teria for a) simulation and b) minimization. Based on earlier discus- 
sions, these criteria are a) N\v^„= N{vp\n and b) N {vs\n ~ 0* 
Using these expressions, and Equations 7 and 8, there results: 
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( 9 ) 



suppression 

In Equation 9, the filtering operation N{ } identifies the portion of 
the upwash velocity wave number spectrum (or eddy field) that contri- 
butes to tone noise. This subject has been addressed in the Phase I 
contract report, ref. (33). Models for the contraction and screen ef- 
fects, denoted symbolically by C [] and B [J respectively, are dis- 
cussed in this report. 

A model for the atmospheric contribution to the external velocity is 
also defined in (33). The development of models for the distortion ve- 
locities generated by stand structure, ground planes, etc. and consid- 
eration of boundary layer and angle of attack effects are outside the 
scope of the current contract, but should be considered in future work. 


Equations 9 and 10 can be used to illustrate the ICS design process. In 
both cases, the ICS designer can control the suppression characteris- 
tics of the screen, i.e. in essence he controls the operation ^[] 
and can also control the location of the screen. The design procedure 
to be developed will be to define the screen characteristics and loca- 
tion in a manner that best satisfies Equations 9 or 10, subject to the 
constraints of minimizing ICS self-generated distortions and acoustic 
transmission losses across the screening. 

5.2 Distortion Generation Constraint 


The ICS should be designed to minimize or eliminate distortions inside 
and outside the boundary layer. Two possible design approaches for 
achieving this can be identified. Firstly, design the ICS without any 
significant distortion generator (such as struts and other inhomogenei- 
ties) in the structure. Alternatively, if such distortion generators 
exist then the distortions generated must be suppressed by a downstream 
device that does not have additional distortion generators or by the 
contraction that occurs downstream of the ICS. 
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5.3 Acoustic Transmission Constraint 

Two approaches are available that account for the acoustic transmission 
characteristics of the inflow control screen. Either it can be required 
that the screen have negligible effect on the radiated sound field or 
it can be required that, in the event that there is an effect, that 
effect is to be defined. The second criterion is the less satisfactory 
on two counts. Firstly, the effect of a screen on the radiated field 
is, in general, a function of the source. Thus acoustic data gathered 
using different engines with the same inflow control screen could re- 
quire different corrections and even data from the same engine at dif- 
ferent operating conditions could require different corrections. Se- 
condly, the necessity of having to correct acoustic data for inflow 
control screen effects increases the complexity of, and thus probabil- 
ity of error in the data reduction process. 

An inflow control screen design that has as small effect as possible on 
the sound radiation field is thus the preferred result. 
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6.0 FAN FACE VELOCITY FIELD DEFINITION (PHASE II, TASK A) 

6.1 Introduction 

During the JNRP, (which was aimed at reducing the noise of a P&WA JT9D 
mounted on a Boeing 747), measurements were taken with pressure trans- 
ducers surface mounted on the rotor blades. The same instrumented JT9D 
was operated at a test stand, with and without an inflow control 
screen, and in flight. It was the objective of this task to use these 
BMT data and hot film data to specify the turbulence velocity charac- 
teristics - intensity and length scale - at the fan face for each of 
these configurations over a range of engine speeds. These characteris- 
tics would then be used to provide an assessment of the preliminary ICS 
design system. In order to quantify the velocity field, it is necessary 
to define the response function relating the pressure sensed by a BMT 
to the unsteady upwash velocity field at the fan face. This response 
function was obtained by measuring the velocity fluctuations in the 
vicinity of the fan face with a hot film, and, by sampling of the BMT 
signal at each revolution, determining the corresponding pressure fluc- 
tuations at the blade surface. However, the once per revolution rate at 
which the BMT signal is sampled imposes an upper frequency limit on the 
empirical response function and so, beyond this frequency, a Sears re- 
sponse form is assumed. This response function can then be applied to 
the BMT data for the different conditions and speeds to determine a 
velocity spectrum, provided the BMT is responding to the velocity 
field. However, there is some question, at this time, as to whether 
the BMT signal is due to the inflow velocity field for all configura- 
tions and speeds. The possibility is discussed that, with the ICS and 
in flight, the BMT's may be responding to other stimuli such as the 
local noise field. Consequently application of the response function to 
the BMT data has been postponed pending resolution of this question. 

6.2 Discussion of Blade Mounted Pressure Transducer Measurements 

The BMT's with which the JT9D rotor was instrumented were thin film 
pressure sensors manufactured by Kulite Corporation. The frequency re- 
sponse of this sensor is greater than 20,000 Hz. The diameter of the 
sensor is .203 cms. while that of its diaphragm is .102 cms. The total 
sensor installation was less than .046 cms. above the blade surface, 
including a thin coating of RTV to provide protection from debris. A 
typical installation is shown in Figure 3. The transducers were located 
on the pressure surface approximately .318 cms. from the leading edge 
of the fan blade. Some were mounted in the tip region of the blade such 
that they were imnersed in the nacelle boundary layer. Others were in 
the free stream region of the fan face. The leadwork was run along the 
pressure surface of the blade and attached to a telemetry transmitter 
in the nose cone of the engine. The data signals were transmitted to an 
antenna on the nacelle wall and subsequently recorded on magnetic tape. 
In order to relate this signal to specific circumferential locations. 
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an impulse was generated every time a particular blade passed a pick up 
mounted on the fan case. This signal was recorded simultaneously with 
the data signal. The BMT data has been analyzed, in part, by digital 
techniques to allow detailed characteristics of the signal to be eval- 
uated. 

The data were digitized such that the signal from each transducer was 
sampled at the rate of 512+1 times during each revolution of the fan. 

At a typical approach condition of 2400 RPM, this corresponds to a 
digitizing rate of about 20,000 samples per secbnd. Data record lengths 
were such that they contained about 1000 revolutions of information. 

The digitized data were stored on tape and input to an analysis program 
that was based on procedures and algorithms developed by Hanson (33). 

This analysis program is capable of separating the deterministic, once 
per revolution part of the BMT signal from the total signal. This de- 
terministic part of the signal is obtained by averaging the value of 
the pressure at a particular circumferential location over the total 
number of revolutions on the record. This is performed at each circum- 
ferential location at which the data was digitized and so the determin- 
istic pressure pattern, defined at 512 angular locations, was con- 
structed. Examples of the average pattern are shown in Figure 4 for the 
three configurations (with ICS, without ICS, inflight at 610 m.) at 
around 2400 RPM rotor speed. The signal from the transducer located on 
the blade, away from the nacelle boundary layer is used. The patterns 
are characterized by an oscillation over one revolution which is be- 
lieved to be due to the presence of a downstream pylon. Superimposed on 
this is a high frequency oscillation due to downstream blade rows. In 
addition, the configuration without the ICS, shows other disruptions 
between 18QO and 30° measured from the top of the inlet counter- 
clockwise looking aft (OCCA) . The variance of the BMT signal about 
this nacelle fixed pattern has also been determined for these cases and 
is shown in Figure 5. This quantity is a measure of the unsteadiness of 
the pressure on the blade at each circumferential location in the na- 
celle fixed coordinate system. The configuration without the ICS in 
place shows greater unsteadiness than either the ICS case or the in- 
flight case, although there is a region from IQO to 120^ CCA where 
all three levels are comparable. When Figures 4 and 5 are compared to a 
photograph of the test stand (Figure 6), it is noted that the region of 
high unsteadiness observed without an ICS in place corresponds in an- 
gular position, to the location of the test stand support structure and 
ground plane. The low unsteadiness region occurs where there is the 
lowest obstruction to the inflow. From these figures several inferences 
may be drawn about the pressure field sensed by the transducers at 
107.75 cm. radius. 

0 Without an ICS in place the unsteady pressure field is domin- 
ated between 120° and lOO CCA by stand/ ground generated 
distortions. 
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Figure 5 The variance of the BUT pressure signal. Engine speed: 
rnm. BUT radius: 107.95 cm. 



0 Without an ICS in place the deterministic pressure field con- 
tains contributions from stand/ground generated distortions in 
the region between 180° and 30° CCA. 

0 With the ICS in place the stand/ground generated unsteadiness 
is substantially reduced. 

0 With the ICS in place, and in flight at 610 m., the unsteadi- 
ness level is approximately uniform circumferentially. 

0 The level of unsteadiness observed with the ICS in place is 
approximately the same as that observed in the low unsteadi- 
ness region with no ICS mounted. 

0 The level of unsteadiness observed inflight at 610 m. is about 
the same as that observed in static operation using the ICS. 

While it is apparent that the high unsteadiness region, observed with 
no ICS, is due to the presence of the ground and stand structure, it is 
not clear what the cause of the unsteadiness is in the low region. It 
is possible that this region is due to the ingestion of atmospheric 
turbulence directly, however, it is also possible that some other un- 
steadiness floor is encountered here; e.g., the local acoustic field. 
This same uncertainty exists also as regards the source of the un- 
steadiness observed in flight and with the ICS in place. 

The ICS design system aspects to be assessed using these data are the 
atmospheric turbulence, contraction and screening models. This assess- 
ment can be performed if the BMT is responding to the directly ingested 
atmospheric turbulence field. However, it is not known at this time if 
this is the case. Consequently, further investigation is necessary to 
resolve this problem. The response function (relating blade surface 
pressure to upwash velocity) required for this assessment has been de- 
veloped and is discussed in the following section. 

6.3 The BMT Response Function 

The fluctuating pressure on the blade surface, in the absence of other 
effects, is due to the fluctuating angle of attack of the flow incident 
to the blade. The nature of this angle of attack variation is deter- 
mined by the unsteady upwash component (that component normal to the 
mean velocity vector) of the incident velocity field. The BMT response 
function, 7(a)) is the ratio of the power spectral densities of the BMT 
pressure signal, P((o) , and the upwash velocity, (m) , producing it; 
i.e.. 


T(.) 


P(ui) 

FnC^) 


( 11 ) 
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The power spectral density of the upwash velocity signal is found from 
the hot film measurements made during the JNRP. The hot film measure- 
ments, made with no ICS in place, were taken in a region where the 
BMT's were responding to the inflow unsteadiness generated by the 
ground and/or stand structure and so an estimate of the BMT response 
function using these measurements can be made. 

The hot films were used to measure the streamwise and circumferential 
components of the unsteady inflow field at a station a short distance 
upstream of the fan. The films were mounted on a probe in the inlet. 

The installation is shown in Figure 7. Some films lay within the na- 
celle boundary layer and others were in the free stream. From the ve- 
locity measurements, the power spectral densities of the streamwise and 
circumferential velocity components, (w) and (w) respectively, 
were found. Using the relationships between the various velocity com- 
ponents, shown in Figure 8, the expression for the power spectral den- 
sity of the upwash velocity component, /^(w) , was determined in terms 
of the power spectral densities, /yCw) and and the cross spectral 

density of the streamwise and circumferential velocity com- 

ponents. The derivation of this expression. Equation 12, is given in 
Appendix A. 


(u) - F, C<j) cos ^ y t Fz (cj) sin^ !^(Ftz (.u)))sinZi 02) 


For the computation of the upwash velocity spectrum, it has been as- 
sumed that the cross spectral term is negligible. Support for this con- 
tention is provided by a typical coherence plot shown in Figure 9 and 
an order of magnitude analysis described in Appendix A. The coherence 
is defined as , 

(w)l 

F/ 

With this assumption Equation (12) becomes 

Fn(i^) = F^ M cos^ y t Fz iu)) Sin^ ^ (13) 

The corresponding blade pressure power spectral density was obtained 
from BMT's located at roughly the same radius as the hot films. Since 
it was not possible to use simultaneous hot film and BMT measurements 
because of the contamination of the BMT data by the wake of the hot 
film probe, the BMT measurements were taken with this probe removed. 

The ambient conditions under which the two sets (hot film and BMT) of 
measurements were taken were different. However, data taken from the 
same BMT under different ambient conditions, without an ICS, indicate 
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that the inflow fields are substantially the same. This data is shown 
in Figure 10. The plots are of the circumferential variation in the 
pressure signal variance (a measure of the unsteadiness of the inflow) 
at two different times. Thus, it may be assumed that the film data was 
gathered under substantially the same inflow conditions as the BMT 
data, even though the measurements were not taken simultaneously. 

Now the BMT pressure data may be viewed in the B plane, where the 
pressure at the BMT is a function of the circumferential location and 
time. As shown in Figure 11 the BMT measures the pressure in this plane 
along lines determined by the velocity of the BMT relative to the mean 
inflow velocity. The dimension is related to an axial coordinate fixed 
with respect to the inflow by the mean inflow velocity. As can be seen 
from Figure 11, only once per revolution does the angular location of 
the BMT correspond to the angular location of the hot film. Further, it 
is necessary to define the temporal behavior of the pressure signal at 
that angular location in order to compute the BMT power spectral den- 
sity required for Equation 11. Therefore, to determine the temporal be- 
havior of the BMT signal due to the flow that passed in the vicinity of 
the hot film probe, the pressure signal must be sampled at this angle 
(202° CCA), during each revolution (Figure 12). At an approach fan 
speed of 2400 RPM, the sampling rate is about 40 samples per second, 
which provides a data point every .025 seconds. Consequently, a reli- 
able estimate of the pressure power spectral density is available only 
to about 20 Hz. To obtain an estimate of the pressure power spectral 
density beyond 20 Hz it was assumed that the blade had a Sears response. 


This assumption leads to the following form for the response function. 


T (lj) 


f\j) 

Frku) 


where ^(co) is Sears' 


/< 20 Ha 

f> 20 Ha 

function. 


(14) 


The response function computation has been performed for two BMT's lo- 
cated at 107.95 cm. radius on two blades. The engine condition was 
about 2400 RPM fan speed and with no ICS. The response functions are 
shown in Figure 13. At the condition for which they were calculated 
there were strong indications that the BMT's were responding to the 
inflow field, so that the calculations are valid. 

The reason for the difference in the two curves is unknown. Conceptual- 
ly, if both blades, transducers and mountings are geometrically simi- 
lar, both should give the same response function. Differences may be 
due to differences in the way the transducers are mounted or to slight 
geometric differences between them. 
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Figure 10 The variance of the BUT pressure signal at different times. 
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The application of these response functions to the BMT data obtained in 
flight and with the ICS in place has been deferred until the questions 
discussed earlier, concerned with determining whether the BMT signal is 
due to inflow velocity perturbations or some other cause, have been 
resolved. 
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7.0 THE EFFECT OF A FLOW CONTRACTION ON CONVERTED TURBULENCE 
(Phase II, Task C) 

7.1 Introduction 

The flow field into an engine is neither uniform nor steady. The un- 
steadiness may be considered to have two parts: (i) that part which 
relies upon the engine flow field for its existence and(ii) that part 
which does not. Included in the first category are such disturbances as 
ground or other surface associated vortices, engine induced wakes and 
the nacelle boundary layer. The second category consists of distur- 
bances of the atmosphere due to the shear of the earth's boundary layer 
and buoyancy forces. These disturbances are referred to as atmospheric 
turbulence and its characteristics depend on such factors as terrain, 
altitude and temperature distribution. It is not isotropic. However, it 
tends to become more so as altitude increases, i.e., more of the tur- 
bulence spectrum has an isotropic form. The atmospheric turbulence phe- 
nomenon and representations of it are discussed in detail in the Phase 
I Report of this contract (34). The effect of the flow contraction on 
this atmospheric turbulence as it is convected into the engine is dis- 
cussed in this section. 

7.2 Background 

The effect of a contracting flow on a turbulent field became of inter- 
est when it was observed that unsteadiness in wind tunnel working sec- 
tions was significantly reduced by an upstream contraction. Prandtl 
(16) analyzed the phenomenon using very simple models that nevertheless 
provided a description, in terms of velocity changes, that corresponded 
well, qualitatively, with subsequent more complicated models. He showed 
that transverse velocity perturbations are amplified by the contraction 
while the streamwise component is attenuated. Taylor (34) in a more 
soundly based analysis showed that the problem could be linearized pro- 
vided that turbulent inertial and viscous forces could be neglected. 

The consequence of his linearization is that the relationship between 
the pre and post contraction vorticity is determined by the change in 
the relative position of the fluid particles and is thus given by 
Cauchy's equation (36). Taylor's linearization constrains the rapidity 
of the flow contraction for a given initial turbulence field. Taylor 
solved the problem for particular convected velocity fields. Subse- 
quently Ribner and Tucker (37) and Batchelor and Proudman (38) using 
Taylor's linearization and a Fourier representation of the turbulent 
velocity field solved the problem for arbitrary homogeneous convected 
velocity fields and contractions. Both sets of authors investigated 
isotropic initial turbulence convected through an axi symmetric contrac- 
tion. Batchelor and Proudman showed that the mean square velocity rela- 
tionships were independent of the initial isotropic spectrum. Ribner 
and Tucker determined expressions for certain post contraction one di- 
mensional spectra, having assumed that the initial form was of the 
Dryden type. Pearson (39) linearized the problem by retaining the vis- 
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cous decay term, and neglecting the turbulence inertial effects. He 
examines, in particular, the effect of an axi symmetric contraction on 
weak isotropic turbulence and shows that as the contraction ratio tends 
to infinity the suppression of the streamwise pertubation velocities is 
accelerated and the amplification of the transverse components is fi- 
nally balanced by viscous decay. In general, however, the total turbu- 
lent energy for a uniform irrotational distortion is found to increase 
without limit. These two linearizations of the Navier Stokes equations 
neglect turbulent inertial forces, with the result that the "return to 
isotropy" tendency observed experimentally is not modelled. 

The attempt by Ribner and Tucker to provide experimental verification 
of the theory by using previously acquired data was not successful, due 
in part, to questions of data validity, and, in part, to violation of 
linearization assumptions. Subsequent comparisons by Uberoi and Hussein 
and Ramjee (40,41) also floundered on this point. In addition, there 
was some possibility of acoustic contamination in Uberoi's data. How- 
ever, both concluded that the theory was most in agreement with_ 
measurement at low contraction ratios. They observed that the linear 
theory, beyond a certain contraction ratio, over predicts both the in- 
crease in transverse velocity perturbations and the attenuation of the 
streamwise component. This observation suggests that turbulence iner- 
tial effects are important in this region. It is significant, as Ribner 
and Tucker note, that the measurements themselves are inconsistent, 
indicating the dependence of the phenomenon on turbulence and con- 
traction characteristics and also the difficulty of making uncontamin- 
ated measurements. 

Townsend (42) also made detailed measurements of turbulence convected 
through a small contraction (4:1) and concluded that, for small 
strains, the structure of the energy containing "eddies" is similar to 
that produced by instantaneous strain. He further detected, at larger 
contraction ratios, an equilibrum condition at which the turbulence 
mean square velocity components are equal, though still changing. 

Tucker and Reynolds (43) disagreed that this equilibrium condition 
existed at the strain that Townsend had determined and further main- 
tained that the total turbulent kinetic energy was reasonably predicted 
by the linear theory if it was corrected for viscous decay. However, 
they admitted that individual component changes were not satisfactorily 
predicted. Both sets of measurements indicated an initially rapid re- 
turn to isotropy after the distortion has ceased. 

The measurements to date have been made in flow contractions of 
velocity ratio that is small compared to that encountered in static 
engine operation. The details of these flows have not been well 
predicted by the linear theory which is supposedly valid only for the 
physically unrealizable instantaneous contraction. 

Based on this review, it was concluded that there is no contraction 
model that has been developed that can be directly applied to the de- 
sign of inflow control structure. 
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It is in part the purpose of the contract to develop such a model by 
building on existing analytical theories and providing empirical cor- 
rections to these theories. The empirical corrections will be derived 
by obtaining data on the effects of contraction on turbulence and mean 
flow distortion characteristics for contraction ratios and contraction 
rates that are more representative of those that result during static 
testing of engines. Since the effects of contraction will be based on 
existing models, a discussion of their underlying theory will be con- 
sidered in the following sections. 

7.3 Theory 

The Navier-Stokes equations provide a suitable starting point for the 
discussion of the problem. Lamb (36) has shown that they may be written 
in terms of vorticity of the flow for an incompressible, constant vis- 
cosity fluid with no body forces. 

^ f y ( 15 ) 

DJ: 

If it is assumed furthermore, that the mean flow is irrotational , then 
the vorticity vector is solely due to the turbulent velocity fluctua- 
tions. Thus after expanding slightly. Equation 15 becomes 

= u, -f u>J t ( 16 ) 

and the dynamics of the vorticity variation may be ascertained. Equa- 
tion 16 states that the rate of change of the vorticity of a point mov- 
ing in a fluid is determined by the sum of (i) the production of vor- 
ticity by the mean motion straining the turbulent "eddies", (ii) the 
transfer of vorticity amongst turbulent "eddies" by one "eddy" strain- 
ing another and (iii) the diffusion of vorticity by viscosity. The 
rapid distortion assumption allows the two latter contributions to be 
ignored on the grounds that the flow is distorted so quickly that no 
vorticity is transferred or dissipated in the contraction time. It is 
this assumption that was invoked by Taylor, et al , to simplify and 
solve Equation 16. 

The validity of the rapid distortion assumption must be assessed rela- 
tive to the flow contractions experienced during engine operation both 
statically and in flight. An order of magnitude analysis of the terms 
in Equation 16, and following Ribner and Tucker is described in Appen- 
dix B. Application of this analysis to the engine contraction problem 
shows that under flight conditions the rapid distortion assumption is 
valid, whereas under static conditions, it is not. 
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Accordingly, use of results based on the simplified form of Equation 16 
should be valid in flight, but questionable for static applications. As 
discussed above, the linearized theories will be used as the basis for 
our models, with empirical corrections used to account for the errors 
introduced by the neglect of the last two terms in Equation 16. 

This simplified equation was solved in detail by both Ribner and Tucker 
and Batchelor and Proudman for an isotropic turbulence field convected 
through an axi symmetric contraction. Ribner and Tucker's treatment is 
discussed in the Phase 1 report. Batchelor's somewhat more elegant 
formulation produced the same results, namely that 
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The length scales shown above are derived in Appendix C. The subscripts 
refer to the vel ocity component and direction of the dimension. For 
instance, is the length scale of the axial velocity component in 

one of the transverse directions. There are nine of such length scales 
in general associated with the three Cartesian coordinates and veloci- 
ties. The above length scale ratios are true, within the assumptions of 
the linear theory, for any initial isotropic turbulence field. The 
length scale ratios not formulated above depend on the particular form 
of the spectrum tensor of the incoming turbulence and therefore cannot 
be generalized. 
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7.4 Eddy Model 

Prandtl used simple "eddy" models to solve the contraction problem. He 
considered two types of flow elements; a steady axial velocity deficit, 
and an axially oriented vortex filament and their modification by an 
axisymmetric flow contraction. Considering the first distortion ele- 
ment, with reference to Figure 14 and applying Bernoulli's incompressi- 
ble equation to two streamlines and assuming that before and after the 
contraction there are no transverse static pressure gradients, yields: 

[/^/9 ” 2 
Now if 
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Figure 14 Contraction of a small steady axial velocity distortion. 
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Thus Prandtl determined that small axial distortions were attenuated 
according to the inverse of the flow contraction ratio. For the convec- 
tion of the axially oriented vortex, again the convecting flow was 
assumed to be axial before and after the contraction. Conserving circu- 
lation around a circular path in a plane perpendicular to the convect- 
ing flow, (see Figure 15), yields 

^CB 

So that (19) 

^( 9 /? r ~ ^ 96 
or 

^9B _ 

U 

Thus from this type of distortion element Prandtl deduced that trans- 
verse velocity components were amplified by the flow contraction ac- 
cording to the square root of the contraction ratio. Now if the turbu- 
lence field is considered to be made up of many eddies and if, further- 
more, the unsteady velocities associated with this eddy field behave, 
through the contraction, in the same manner as a velocity field made up 
of many of Prandtl 's elements, then an ensemble average over the up- 
stream and downstream velocity fields, yields relationships between the 
mean square values of the perturbation velocity components. 




J_ (tfjL 4'J)/. _ J 

I," ’ ' 


Furthermore an examination of the distortion of a cubical fluid element 
by an axi symmetric contraction where the field particle displacements 
are dominated by the mean flow gradients (Figure 16) produces estimates 
of the eddy length scale changes. 

L-m Lsa ( 21 ) 


These length scale ratios are considerably simpler than those derived 
from Taylor and Batchelor's work. Here, there is no dependence on the 
velocity component and the ratios are only functions of the direction 
denoted by the subscript. The length scale ratios then for all three 
velocity components in a given direction are the same. 
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Figure 15 Contraction of a steady axial vortex. 
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Figure 16 The axisgimetric contraction of a fluid element. 
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Expressions 20 and 21 describe the effect of an axi symmetric contrac- 
tion on a turbulent field. In view of the absence of viscous decay ef- 
fects on the distortions elements, together with the assumed absence of 
any interaction between these elements, this model can be considered to 
be an approximation of the general solution of Taylor's linearized 
problem, determined in Refs. 37, 38. A comparison of the two sets of 
results is shown in Figure 17. As can be seen the trends are similar 
for the mean square values of the velocity components, and for the 
larger contraction ratios. The corresponding predictions of the kinetic 
energy ratio across contraction are shown in Figure 18. The comparison 
of scale changes predicted by the two theories is shown in Figure 19. 
The dashed lines are the scale ratios given in Equation 21 while the 
solid lines are the ratios obtained from Batchelor and Proudman's re- 
sults. 

The predictions of the scale change in Ln are not in close agreement 
and this reflects the different upstream turbulence fields that were 
examined. The general significance of the integral length scales is 
somewhat obscure on a physical basis, they are however, a measure of 
the one dimensional spectral shapes. The one dimensional power spectral 
density Fz is the power spectral density of the ith velocity 

component as a t unction of the wavenumber component in the jth direc- 
tion. This power spectral density is analogous to the Fourier decompo- 
sition of a particular velocity component in a periodic velocity field 
in a particular direction. 

If in addition to the earlier assumptions it is required that the 
"eddies" be axially oriented, then the downstream one dimensional spec- 
tra are functions only of the corresponding upstream one dimensional 
spectra and the contraction characteristics. In this situation the en- 
ergy in the wave number range about will be converted into 

energy in the range A -^^3 , see Figure 20. Now with the one 
dimensional spectra represented by 
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Figure 19 The effect of a contraction on the turbulent scale ratios. 
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Where uj_- u,, are the Cartesian components of the un- 
steady velocity field and are the Cartesian 

wave number components. 

Comparisons between two downstream spectra as predicted by Equation 22 
and Ribner and Tucker are shown in Figures 21 and 22. In view of the 
discrepancy noted previously for the Ln integral length scale, the 
large difference in the corresponding spectra is not unexpected. How- 
ever, there is close correspondence between the predictions of the 
transverse spectrum, again in accordance with the length scale compari- 
sons. 

The predictions using the eddy model are then close enough to those 
using the more rigorous theory to provide a basis for experimentation. 
The simplicity of the one dimensional spectral relationships recommends 
their use in the inflow control system design procedure. Figure 23 con- 
tains a summary of the eddy model results. 

The linear theories account only for the distortion process and this 
proceeds, according to Uberol , et al , unrelieved by an turbulence 
inertial or viscous decay effects. Consequently the higher the 
contraction ratio, the more in error will be the linear theory. 

7.5 Conclusions 

I. Based on the discussions in the section, the simple theory of 
Prandtl , summarized in Figure 23, will be used as the basis for 
the models that account for contraction effects on turbulence. 
Tendencies toward isotropy, observed experimentally for small con- 
traction ratios relative to the static case, will be incorporated 
empirically as necessary, by use of the test data from Task F of 
this contract. 

II. The measurements of Uberoi , et al , indicate that the contraction 
process is characterized by: 

0 A distortion of the turbulence field. 

0 Amplification of the transverse turbulent velocity components . 

0 Attenuation of the streamwise turbulent velocity components, 
at least, initially. 

0 An equilibrium condition attained after a certain amount of 
contraction beyond which further distortion does not occur. 
This equilibrium condition is determined by the contraction 
rate variation through the contraction. 

0 A rapid initial rate of return to isotropy followed by an 
asymptotic approach to that state. The initial return rate is 
determined by the distortion of the field. 
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Figure 23 The eddy model for the contraction of turbulence. 
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III. The comparison of the rigorous and simple theories should be 
viewed with this in mind. The main points of the comparison are: 

0 The predictions of the turbulent kinetic energy ratio across 
the contraction produced by both theories are in close agree- 
ment. 

0 The predictions of the ratio of the RMS values of the indi- 
vidual transverse velocity components across the contraction 
produced by the two theories are in reasonably close agreement 

0 The predictions of the ratio of the RMS values of the stream- 
wise velocity component across the contraction produced by the 
two theories are not in good agreement, although the trends 
with increasing contraction ratio are similar. 

0 Correspondingly the predictions of integral length scales and 
one dimensional spectra of the transverse components of turbu- 
lent velocity produced by the two theories are in reasonably 
close agreement. But the predictions for these quantities for 
the streamwise component of turbulent velocity are not in good 
agreement except at low contraction ratios. 

In general then, the more detailed the description of the phenome- 
non, the more divergence between the predictions of the simple and 
rigorous theory. 

IV. The applicability of the linear theories to a JT 9 D operating in 
flight and statically was examined and it was concluded that: 

0 The linear theory describes the inflight contraction of atmos- 
pheric turbulence for the JT 9 D. 

0 The linear theory does not describe the contraction of atmos- 
pheric turbulence by the JT 9 D at the test stand. In this situ- 
ation the non linear turbulence inertial effects are not neg- 
ligible and the linear theory will overpredict the distortion 
of the turbulence field. 

V. The contraction ratio in flight is of the order one and statically 
this ratio is large and so, bearing the comparison of the linear 
theories in mind, it is suggested that: 

0 The simple model is quite adequate to describe the contraction 
of turbulence in flight for the JT 9 D. 

0 Either the simple model or more rigorous theory produce an 
upper bound to the distortion of the atmospheric turbulence 
produced by the contraction due to the JT 9 D operating sta- 
tically. A lower bound is provided by the assumption of a com- 
plete return to isotropy. 
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8.0 THE EFFECT OF A SCREEN ON ATMOSPHERIC TURBULENCE 
(Phase II, Task C) 

8.1 Introduction 

The problem of the interaction of screens and flow was originally ad- 
dressed to optimize the design of screens for use in wind tunnels. 

These screens had been observed to suppress turbulent fluctuations in 
wind tunnel flow. The solution was provided, after several partially 
successful attempts by other investigators (16, 17, 18), by Taylor and 
Batchelor (19). They solved the problem for thin gauzes subject to 
flows of small turbulent intensity. They deduced the relationships be- 
tween the upstream and downstream intensities and Batchelor (43), in a 
later, more elegant analysis showed the relationships between the 
velocity spectrum tensors of the two flow fields. They also considered 
specifically the case of initially isotropic turbulence and showed for 
this case, that the intensity ratios are independent of the form of the 
energy spectrum. The analysis, as in the contraction problem, was made 
tractable by a linearization of the Navier Stokes equations. This solu- 
tion, as in the case of the contraction problem, results in no repre- 
sentation of the inertial "return to isotropy" phenomenon or viscous 
decay. 

As noted before, the ability of the screen to attenuate unsteadiness in 
the flow is determined by two of its characteristics; (i) the pressure 
loss coefficient and (ii) the amount of turning the screen inflicts on 
the flow. The screens remove vorti city from the flow by absorbing it 
as a torque. In the following sections, a brief survey of the theory of 
Taylor and Batchelor is presented together with their approximate form. 
Then, using Lumleys results, it is shown that the results of Taylor and 
Batchelor's analysis for gauzes are applicable to honeycomb screens, if 
the incident turbulence scale is large compared to the honeycomb 
length. Subsequently, expressions for determining the pressure loss 
coefficient, /C , and flow angle ratio, a, for honeycomb and perforated 
plate are given. 

8.2 General Theory for Thin Screens 

The fundamental assumption of Taylor and Batchelor's approach is that 
the turbulent velocity fields upstream and downstream of the screen do 
not decay and the distribution of energy is unchanged away from the 
screen. This, (coupled with the knowledge that the effect of the screen 
on the adjacent flow is due entirely to the pressure change across it, 
and is therefore irrotational ), allows the velocity field induced by 
the screen to be superimposed on the turbulence field. In short, the 
convection process is linear, the effect of the screen on the flow 
field is linear, and therefore, the velocity fields can be summed. The 
flow is considered inviscid, except insofar as the screen constitutes 
an impedance to it, thus the turbulence generation characteristics of a 
screen are not accounted for. Furthermore, in Taylor and Batchelor's 
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treatment, the screen afflicts the flow with a sudden pressure drop and 
a sudden change in direction. The screen is therefore considered in- 
finitely thin and appears to the turbulence as a continuum requiring 
that the transverse turbulent scales are large compared to the mesh 
size. 

The condition under which the decay processes are negligible is that 
the time scale of the decay is large compared to the time taken for a 
particle to traverse the screen influenced flow region so that 


J— 


a 


u 



( 23 ) 


As in the contraction problem Batchelor approximates the derivative 
with an empirically derived relationship 

Cbt t ( 2 ^) 


to reduce the criterion to 





( 25 ) 


Again, this empirical expression may not be valid for the screen region 
in view of the disruption of the decay process there. It does, however 
provide a tractable criterion for linearization and, in the absence of 
any superior relationship, will be used here. It should be noted that 
this linearization, performed in the absence of mean flow gradients, 
leads to turbulent velocity fields where vorticity is conserved. Only 
in the vicinity of the screen itself is there a change in vorticity. 

Taylor and Batchelor argue that the effect of the screen on the flow is 
irrotational and therefore the contribution to the velocity field due 
to the screen may be expressed as a velocity potential. Thus the form 
of the total unsteady velocity field is 


U(X,- U,t, X2,X^) -t (Xj£) ( 26 ) 


55 



where u is the turbulent velocity vector far away from the screen. 

With this velocity field form, and the linearization discussed pre- 
viously, together with the small perturbation assumption, the problem 
is solvable since the upstream and downstream velocity fields may be 
related at the gauze through the screen pressure loss coefficient. 



The remaining relationships are those produced by consideri 
0 The Turning of the Screen 



0 Continuity Across the Screen 



0 Continuity in the Upstream and Downstream Flow Fields 
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Local outflow angle 
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when the inequality (25) is satisfied, i.e., for small perturbations. 
The coordinate system is shown in Figure 24. 

If the turbulent velocities arid potential are Fourier decomposed, then 
this set of equations allows the statistical characteristics of the 
upstream and downstream velocity fields to be determined in terms of 
the pressure loss coefficient,/^', and the flow angle ratio a . 

For the particular case of initially isotropic turbulence, Taylor and 
Batchelor found that the ratios of the mean square values of the up- 
stream and downstream velocity components away from the screen were the 
same for all isotropic fields. 

Specifically they found that 
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It is also possible to show (see Appendix D) that the form of certain 
integral length scale ratios is independent of the initial turbulence 
spectrum. Specifically 
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For initially isotropic turbulence these ratios may be evaluated using 
Taylor and Batchelor's expressions for mi and|^ for a gauze. For 
this case a and K are apparently related so that the velocity ratios 
may be plotted against K as shown in Figure 25. The corresponding 
scale ratios are shown in Figure 26. As Taylor and Batchelor note, in 
general, there exists a screen such that the integral length 
scale Lfis is totally suppressed. This occurs when 


oc/<'-/ -OC =0 


(35) 
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Figure 24 The coordinate system for Taylor and Batchelor's theory. 


58 





The other scales in Equation 34 are suppressed when ot = 0. 

The results of Taylor and Batchelor's analysis have been compared to 
data gathered by Townsend (45). On a kinetic energy basis there is 
close agreement between measurements and predictions (Figure 27) 

(Ignore the square symbols at this time). Comparisons of the mean 
square velocity ratios highlighted the absence of any "tendency to 
isotropy" considerations in the theory, Figure 28. (Ignore dashed lines 
at this time.) Townsend's measurements indicated a less anisotropic 
downstream field than Taylor and Batchelor would predict. 

It is possible to determine, if a form for the upstream isotropic spec- 
trum tensor is assumed, the changes in the one dimensional power spec- 
tral densities as a result of the turbulence being convected through 
the screen. The manipulations required are complicated and in view of 
the return to isotropy observed by Townsend would produce results that 
require empirical correction. Consequently a simpler model will be de- 
scribed that readily permits calculation of the change in the one di- 
mensional spectra. While these results will probably also require em- 
pirical correction, it is considered that the amount will be no more 
than for the results of Taylor and Batchelor's more rigorous analysis. 

8.3 Simple Screen Model 

In a manner analogous to Prandtl's treatment of the contraction phe- 
nomenon, Taylor and Batchelor developed a simple model that considered 
the effect of a screen on simple flow distortion elements. These ele- 
ments are (i) a steady distortion of the streamwise velocity component 
that varies sinusoidally in the transverse direction and (ii) an arbi- 
trary steady transverse velocity distortion. 

To determine the solutions of the problems, the Navier Stokes equations 
are linearized by requiring that the amplitudes of the distortions be 
small compared to the mean flow velocity. Viscous effects are also ne- 
glected except insofar as they give rise to the resistance of the 
screen. The assumption of incompressibility, applied to the first dis- 
tortion element, leads to the equations of motion of the fluid, up- 
stream and downstream of the screen, in terms of the stream function; 


f— f 2 

V ^ sinkx^. (36) 

ipg = ^ sink. >>2. 


61 



KINETIC ENERGY RATIO 


1.0 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 



6 

o 



A TAYLOR AND BATCHELOR 
O SIMPLE THEORY 
O TOWNSEND'S MEASUREMENT 



K (LOSS FACTOR) 


Figure 27 The comparison between measured and predicted turbulence 
kinetic energy ratios across a gauze. 
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Figure 29 depicts the geometry of the problem. If it is further assumed 
that the fluid perceives the screen as a continuous body then the solu- 
tions of these equations together with the boundary conditions at the 
screen, Equations 38, 39, and 40, allow the effect of the screen to be 
determined. 

0 Turning of the Screen 


[li-L.] = ocj-iti-j 

l<^*, <x, = .o (38) 


0 Continuity Across the Screen 


(39) 

<3^2 = I ’x,=i -0 


0 Pressure Drop Across the Screen 



(40) 


The boundary condition. Equation 40, also makes use of the small per- 
turbation assumption. Solving these equations simultaneously Taylor and 
Batchelor concluded that 


^/a _ /- a-a K 

" l+ct-hK 


From this it can be deduced, in view of the linearity of the equations, 
that any combination of sinusoidal streamwise perturbation will be 
modified according to Equation 41. 

The effect of the screen on this type of distortion is such that in the 
vicinity of the screen, the potential field results in a static pres- 
sure variation in the transverse direction. Fluid is thus directed 
laterally from the high to the low pressure region. Thus a transverse 
velocity component is induced in the flow. The magnitude of this com- 
ponent is determined by the resistance of the screen to the fluid. The 
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flow with the induced lateral velocity component is transported through 
the screen. The screen attenuates the lateral component according to 
the side load applied to the screen. The residual streamwise perturba- 
tion at the screen is attenuated according to the resistance of the 
screen. Downstream as the potential effect of the screen decreases the 
static pressure becomes uniform laterally and the transverse velocity 
component tends to zero. 

The attenuation of a purely streamwise distortion, because of the 
screen induced lateral velocities, becomes then a function not only of 
the screen resistance but also of the ability of the screen to with- 
stand a sideload. 

As can be seen from Equation 41 there exists a condition when the 
streamwise perturbation is totally suppressed, this occurs when 


I i-cx - C< K-= 0 


(42) 


and corresponds to the suppression of the integral length scale Ln 
noted in the previous section. 

In the case of an arbitrary steady distortion in the transverse veloci- 
ty field, the linearization leads to the deduction that the transverse 
velocity field is unaffected by the screen in that no potential field 
is superimposed on the flow by the screen. The change in the flow field 
across the screen is solely determined by continuity and thus by the 
ability of the screen to turn the flow. Thus 


U 


ZB 


u 


= a 


2/9 


(43) 


The resistance of the screen is irrelevant for this type of disturbance. 


Again the linearization of the problem allows Equations 41 and 43 to be 
applied separately to the individual components of a compound velocity 
field and the results recombined to form the downstream velocity field. 
The criteria for this linearization of the Navier Stokes equations are 
roughly 


Uf 


(44) 
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If a turbulence field is considered to be approximated by some super- 
position of the elements discussed above then an ensemble average of 
the square of the upstream and downstream perturbation velocity fields 
produces 



/ /y- a - o A' \ 

I 

(45) 

u}, < 

(46) 
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In order to further develop the simple model it is not necessary to 
evaluate Equations 45 and 46 for the purpose of comparison with predic- 
tions from the more rigorous theory and Townsend's data (45). Figure 27 
shows this threeway comparison for overall turbulent kinetic energy 
change across the screen. In this comparison, a gauze is used for which 
it can be demonstrated that a relationship exists solely between oe 
and K . An implication of a relationship existing between ot and 
for gauzes is that both parameters are determined by the same geometric 
characteristics in the region where the relationship holds. This rela- 
tionship is described and discussed further in Section 8.5. The pres- 
sure loss coefficient, K , is chosen as the independent variable for 
this comparison. The agreement between the two prediction methods and 
the data is generally very good. 


The ratios of the mean square values of the individual velocity com- 
ponents are not as well predicted by either theory, as shown in Figure 
28, although the theories are in reasonable agreement with each other. 
The reason is that downstream of the gauze turbulence inertial effects 
are not negligible and thus there is a tendency for the velocity field 
to become isotropic. Thus a possible modification to the simple model 
would be to assume that the downstream field is isotropic. Such an as- 
sumption, if the incident field were also isotropic would result in the 
relationship 



(47) 


The appropriateness of Equations 45 and 46 or 47 will be examined in 
Section 8.6. 
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The tendency to isotropy also alleviates the distortion of the length 
scales and the corresponding one dimensional spectra. To account for 
this effect, the simple model suggested here will incorporate this ef- 
fect by having the same length scales upstream and downstream and the 
one dimensional spectra will be unchanged in shape, so that 



( 48 ) 




( 49 ) 


Where Aji is determined by Equations 45, 46. 

There is some evidence in Townsend's data that the spectral relation- 
ships are approximately correct. Additionally, at least some of the 
length scales are substantially unchanged by the screen according to 
the more rigorous theory, see Figure 29. 

This simple form of the theory provides a basis for generating empiri- 
cal corrections from the measurements to be taken in the test program 
of Phase II, Task F. It should be emphasized that the theory above does 
not account for generation of turbulence by the screen. This aspect 
will be addressed in Section 8.7. 

The theories described above assumed an infinitely thin screen, and as 
such should be applied only to thin screens and perforated plate. It is 
possible to extend these theories so they are applicable to finite 
depth screens such as honeycomb. This is the subject of the next sec- 
tion. 

8.4 The Application of the Thin Screen Models to Honeycomb Screens 

The analysis of Taylor and Batchelor assumes that the flow experiences 
an instantaneous total pressure drop at the screen. The consequence of 
this model is that the screen is, in effect, considered to be infinite- 
ly thin. Lumley (20) has attacked the problem for a specific situation 
of a honeycomb screen and so a formulation is available for a finite 
thickness screen acting on turbulence. If the contention that the gauze 
theory can be applied to honeycomb is correct, then a limiting case of 
Lumley's actuation should give the same results as the Taylor-Batchelor 
theory. This comparison is made in the following text. 
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Lumley has derived a formulation for the effect of a honeycomb screen 
on a turbulent inflow for a specific set of conditions. He assumes that 
the inflow field is isotropic, that the honeycomb destroys all lateral 
velocity components and finally that downstream of the honeycomb there 
is a return to isotropy. From his theory he computes values of the ra- 
tio of the rms velocities, upstream and downstream, as a function of 
loss factor and the incoming streamwise length scale to cell depth ra- 
tio (Figure 30). 

As this latter quantity becomes larger, the ratio of the mean square 
values of the velocity components tends to an asymptotic form solely 
dependent on the pressure loss coefficient. It is this form that is 
compared with the ratio determined from the simple model, since in this 
limit the honeycomb is perceived by the turbulence as being thin. Thus 
from Equation 46 setting <x=0 to characterize the suppression of the 
transverse velocity by the honeycomb. 


ilM- = / 

uZi 30 +K)^ 


This corresponds to Ti of Lumley' s formulation. A comparison of the 
results from the two theories is shown in Figure 31 and there is close 
correspondence. The asymptotic form of Lumley' s theory is reached for 
streamwise length scales approximately equal to ten times the honeycomb 
thickness. This, therefore, constitutes a possible constraint on the 
application of Equation 50 to honeycombs. 

Having established the region of validity of the simple theory, the 
evaluation of a and K for gauze and honeycombs will now be 
examined. 

8.5 The Pressure Loss Coefficient and Flow Angle Ratio for Thin 
Screens and Honeycomb 

In general, both the pressure loss coefficient, Af , and flow angle ra- 
tio or , are dependent on the geometry of the screen and the nature of 
the flow i .e. 


K - K (Hl , M ,&!,&) (51) 

® = <‘iJU , M , &i,9) (52) 
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Figure 30 The graphical form of the results of Lumleg's theory of 
honeycomb . 
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Where G represents some geometric qualities. If compressibility ef- 
fects are negligible and small incident flow angles are assumed, and if 
further, the turning of the screen is solely due to its potential effect 
on the flow then Equations 51 and 52 become 


k-k (Ri., s;) 

(53) 

a= aC&i) 

(54) 


Now the measurements of Dryden and Schubauer (46) lead to the conclu- 
sion that over the range of their experiment, the phenomenon is only 
weakly dependent on Reyonds Number, and if Equations (53) and (54) are 
examined in this light it seems quite reasonable to suspect that a re- 
lationship exists between a and A' . Taylor and Batchelor determined 
that such a relationship existed from the above mentioned data and 
found that 


« = U f^y , -] ( 55 ) 


was reasonably good fit, Schubauer, Spangenberg and Klebanoff (47) de- 
termined that for K < J their data was well fitted if 



It should be emphasized that these expressions are valid only for the 
small incident flow angles and were determined using gauzes. It is not 
known if other thin screens behave in this manner, however, it will be 
assumed here that they do and that Equations 55 and 56 apply to perfor- 
ated plate screens in particular. 

The pressure loss coefficient, , depends on such geometric parameters 
as the hole density, hole diameter and thickness of the screen. Baines 
and Peterson (48) conducted extensive measurements on various thin 
screens forms including perforated plates and they noted that 


Where Cc is the contraction coefficient of a hole and s is the sceeen 
solidity. Their results were explainable in terms of this expression. 
The contraction coefficient is determined by the detailed geometry of 
the holes, ranging in value from .7 for sharp edges holes (thin perfor- 
ated plate) to unity for well rounded holes. Figure 32 shows how the 
coefficient varies with solidity for different screen types. No depen- 
dence on a Keynolds Number based upon wire diameter or the equivalent 
dimension for a perforated plate was observed between values of 10^ 
to 10^. At low Reynolds Numbers, there is a dependence on that 
parameter, however, it is not anticipated that an inflow control screen 
will operate in this region. 

To find a formulation for a and K that is suitable for a honeycomb 
screen requires consideration of the factors separately. 

For the determination of a , it is noted that irrespective of the 
orientation of a flow with respect to the honeycomb, the flow perceives 
it as approximately a series of flat plates in cascade. Thus, it is 
suggested that the honeycomb acts on the flow in a similar manner to an 
infinite cascade of flat plates at zero stagger angle. The basic mech- 
anism for deflection of the flow, namely, the loading of the material, 
is the same in each case and the magnitude of this loading is deter- 
mined by the pitch-chord or diameter to length ratio. While the cascade 
effect is, essentially, a two dimensional phenomenon and, that due to 
the honeycomb, three dimensional, it will be assumed that the action of 
the honeycomb on the flow is axisymmetric and determined by the two 
dimensional model. 

The solution to the cascade problem is presented by Von Karman and 
Burgers in (49) and is quite simple in form, 

L 

a = d ( 58 ) 


With being the length to diameter ratio of the honeycomb. 

This expression for the flow angle ratio for honeycomb does not have to 
stand alone, for as Prandtl notes, "a honeycomb length to diameter ra- 
tio of two is sufficient to straighten an incident flow field". More 
recently, the finding of Ginder, Kenison and Smith (33) that honeycomb 
of length to diameter ratios greater than two does not produce more 
turbulence attenuation is totally consistent with this form for « , 
since beyond this ratio, a is effectively zero. On a more personal 
level, anyone may blow through honeycomb at quite large angles of inci- 
dence and note that the flow is straightened for quite small values 
of l]cL, 
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Figure 32 Thin screen loss factor. 
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To determine the pressure loss coefficient for honeycomb, Blasius* 
formula for the skin friction coefficient of turbulent flow in a pipe 
will be used, so that 
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This empirical relationship is discussed by Schlichting (50) and is 
accurate in the range of diameter Reynolds Numbers from 3.1Q3 to 
10^. Predictions of honeycomb loss coefficient using this expression 
were compared with measurements made by Paterson (51) (see Appendix E) . 
Agreement was quite close. 

It is quite possible that honeycomb screens will be operating at lower 
Reynolds Numbers than S.IO^ and consequently the laminar pipe resis- 
tance coefficient might provide a more accurate expression. 
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This result of the solution of the Navier Stokes equations for 
Hagen-Poiseuil le flow has compared well with measurements made in pipes 
containing laminar flow. 

It is interesting to compare the effects on turbulence of thin gauze 
type screens with honeycomb. To facilitate this, an expression relat- 
ing OL and K for honeycombs can be obtained by combining Equations 59 
and 60. 


a ss sxp 


‘3IUH 


(61) 


If a Reynolds Number of 1000 is assumed, (which corresponds roughly to 
a 6mm. honeycomb operating in a flow of 2 m/sec), it is immediately 
apparent that a tends rapidly to very small values for increasing K . 
The consequences of this behavior are shown in Figure 33 with predic- 
tions made using the simple theory. The transverse ratio moves rapidly 
to very small values and correspondingly, the streamwise ratio tends 
rapidly to decay as (//■/<)"*. The honeycomb then, in a sense, acts on 
the turbulent field, in an opposite manner to that of the thin screen. 
In the former, the transverse velocity components are preferentially 
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suppressed while in the latter screen type, the streamwise components 
are attenuated most for a given screen loss coefficient. These comple- 
mentary effects raise the interesting possibility of a combination 
screen that can attenuate an isotropic field and produce immediately 
downstream, another field that is at least isotropic as regards the 
mean square values of the component velocities. 

8.6 The Application of the Simple Screen Model to a Screen Operating 
on Atmospheric Turbulence at the Test Stand 

An important assumption that the linear theory rests on is that 



(62) 


This assumption corresponds to the requirement that turbulence inertial 
and viscous decay effects are unimportant whilst the flow is convecting 
through the region of influence of the screen. The validity of this 
assumption can be assessed using the results of the atmospheric turbu- 
lence^odel described in the Phase I Report. A typical value for 

in the atmosphere near the ground is .3, however, at the 
screen, there is a superimposed velocity field due to the engine opera- 
tion. For a JT9D at approach and a screen mounted at 4m. radius, this 
velocity will be about 7 m./sec. The typical mean atmospheric velocity 
used to compute a value of .3 was 2.68 m./sec., so that when the engine 
inflow field is imposed on the atmosphere, an overall intensity on the 
average of .3 x 2.68/4.84 = .167 at the screen would be expected. Con- 
sequently, the linear theory criterion is marginally satisfied for a 
static test situation. In addition, the satisfaction of this criterion 
means that the incident flow angles are sufficiently small so that the 
constancy of a and K is assured and the expressions for these factors 
are valid. The flow does not suffer from compressibility in this re- 
gion. The change in the turbulence field across the screen is thus 
probably reasonably well approximated by the linear models. However, 
the subsequent tendency to isotropy which apparently is greater accord- 
ing to the extent of the anisotropy of the field is not modeled. 

In view of the nature of the return to isotropy as described by Town- 
send (45) in which there is a rapid readjustment followed by a very 
slow return to isotropy, it is probably that neither the anisotropic 
nor the isotropic form of the simple model will exactly describe the 
phenomenon. The anisotropic form of Equations 45 and 46 will over esti- 
mate the distortion of the turbulence field by the screen while the 
isotropic form of Equation 47 will underestimate that distortion. Hence 
the two forms put bounds on the values of the turbulent velocity ra- 
tios, and can be used in this manner. The data gathered in the Phase II 
task of test program will be used to select the model that best des- 
cribes the data and also to provide an empirical correction to that 
model . 
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8.7 Distortions Induced By Screens 

A simple model has been presented that describes the phenomenon of tur- 
bulence convected through gauze type and honeycomb screens. The model 
does not consider self-generated distortions of which there are two 
types: one generated by unsteady deflection of the screen and the other 
consisting of the wakes and vortices shed from screen elements. The 
subject of screen deflection is discussed first. 

If the screen is exposed to an unsteady vortical flow and deflects 
under the corresponding unsteady applied torque, then, as the applied 
torque varies, so will the screen torsional deflection thus imparting a 
swirling motion to the flow. The screen, then, in addition to suppress- 
ing vorticity in the flow, can also add vorticity. This effect can be 
important if the screen deflects significantly as an eddy passes 
through it. For after the eddy has passed, the screen will unwind and 
thereby impart a swirl to the flow that is coherent roughly over the 
area of the initial eddy. Thus eddies having scales that result in sig- 
nificant turbulence rotor interaction noise levels may not be effec- 
tively suppressed. It is, therefore, necessary that this type of screen 
induced distortion be minimized by either making the screen sufficient- 
ly tors ion ally rigid or mounting another screen downstream of it. 

The other type of screen- induced distortion, that due to the wake of 
the screen, in general, will not be significant as regards turbulence - 
rotor interaction noise in view of the very small scales involved. 
However, wakes and vortices shed from screen structural support members 
may be important. This latter distortion can be effectively minimized 
by design or by suppression after it has been generated. 

8.8 Conclusions 

The theory developed by Taylor and Batchelor for the effect of a screen 
on small turbulent fluctuations has been discussed. Their linear theory 
does not account for turbulence inertial effects in the vicinity of the 
screen and thus the tendency to isotropy observed by Townsend and 
others is not described. A simplified form of the results of Taylor and 
Batchelor produces predictions that are in reasonable agreement with 
the more rigorous theory predictions. This simplified model relates 
pre-screen and post screen turbulence characteristics in terms of the 
flow angle ratio a , and pressure drop coefficient K , of the screen 
and these relations are 

0 The turbulence kinetic energy ratio 
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0 The component mean square velocity ratios 




afl ~aK 
l-hatK 




0 The length scale ratios 



B 



0 The power spectral density ratios 

Taylor and Batchelor's theory is applicable for situations where the 
flow perceived the screen as a homogeneous device, i.e., where turbu- 
lence scales are large compared to a characteristic transverse screen 
dimensions. Also a thin screen is assumed. It was demonstrated that the 
model is applicable to honeycomb type screens if they appear thin to 
the turbulent flow field, i.e., the turbulent scales are large compared 
to the honeycomb length. The forms of a and K are different for 
these two types of screens. The expressions identified for these 
quantities are: 


0 For gauze screens 

or perforated plate 
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a= 8-K 

S-K 

K<.1 


where Cc is the hole contraction coefficient and S' is the screen 
solidity. 
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0 For honeycomb screens 


K- 


.3IUH 


A 

cL 


turbulent 



IV A 
Rjl d 



laminar 


where ^ is the length to diameter ratio of the honeycomb. 

The expressions for honeycomb pressure drop coefficient are solely 
based on pipe friction. It is possible that there is a significant con- 
tribution to this quantity from other components of the honeycomb form 
drag. 

The distortions generated by the screen itself were also discussed and 
in summary these are: 

0 Small scale turbulence determined by the detailed nature of 
the screen and insignificant as regards the generation of tone 
noise. 

0 Wakes and vortices generated at screen support structure and 
at screen discontinuities that are possibly significant in the 
generation of tone noise. 

0 Distortions caused by deflection of the screen, upon removal 
of torques applied to the screen by upstream unsteady distor- 
tions. These could be of importance in tone generation. 

The latter two self induced distortion types can be readily avoided or 
eliminated by proper design. 
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9.0 THE EFFECT OF A FLOW CONTRACTION AND SCREENING ON A 
STEADY DISTORTION (Phase II Task D ) 

9.1 Introduction 

Atmospheric turbulence, as discussed previously, is not the sole source 
of distortion of the inflow to an engine. The engine itself induces a 
velocity field over adjacent surfaces and thereby produces distortions 
to the flow such as wakes and vortices. Also the general asymmetry of 
the test stand situation results in a non-uniform inflow to the engine. 
These self-generated distortions are somewhat unsteady, due in part to 
the viscous effects at the generating surface and in part, due to modu- 
lation of the distortion by atmospheric turbulence and engine speed 
variations. Furthermore, the inflow distortion field is inhomogeneous 
and thus when it is modulated, it is perceived by a fixed transducer to 
be non-stati onary. The modulating flow is at low frequency, so that 
over some period of time, the flow over the generating surface may be 
considered to be steady. Therefore, approximately, the effect of the 
modulation is to change the position of the distortion but not its 
structure. Thus in a frame of reference attached to the distortion, the 
velocity field is deterministic, approximately, and the mean value of 
the distortion velocity field exists. The phenomenon so described, may 
be termed "quasi-steady" in that the nature of the distortion as viewed 
from a frame attached to itself is approximately the same as that of a 
distortion from the same generating surface in a steady flow. 

The vortex attached to the ground plane may be considered quasi-steady. 
It does exhibit intermittent growth and destruction (51), however, its 
structure is well defined over periods of time long compared to a rotor 
period and may therefore be treated as a steady distortion. 

This section is concerned with the modification of the mean value of 
the velocity field associated with these types of distortion. The mod- 
els used are substantially those described in Sections 7 and 8. 

9.2 Contraction of a Steady Distortion 

In Section 7, it was determined that distortions of the streamwise ve- 
locity component were attenuated, by a contraction, according to 


Mi« - J- 

AU,. 'A (63) 


The steady velocity component of a wake will in addition be modified by 
dissipation through viscosity over the length of the contraction. The 
attenuation of Equation 63 is therefore a lower bound on the amount by 
which a wake will be suppressed. The contraction of a dissipative wake 
will now be considered. 
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The characteristics of the wake at the fan face of most interest from 
the point of view of fan noise generation are its width and velocity 
deficit. These characteristics are assumed to evolve according to 


A u, <c 



(64) 




which corresponds to the development of the far wake of a rod of diame- 
ter a and drag coefficient . 


Further, to permit a tractable model to be constructed using existing 
theories for wake contraction and dissipation, it is assumed that this 
wake develops through the contraction in alternating purely dissipative 
regions and instantaneous contractions. In this manner, the contraction 
characteristic can be approximated by a series of steps, see Figure 34. 
Each dissipative region is characterized by a value of > Cj^a. ,and 
X . The latter quantities in any dissipative region can be considered 
to be the virtual characteristics that, in a uniform flow, will produce 
the same wake as that present in the dissipative region. In the limit 
of an infinite number of dissipative regions of infinitesimal extent 
Cga and X , like , become continuous functions of position. In 
this limit, the wake characteristics at To can be related to some 
reference wake characteristics at such that 
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( 66 ) 


Now defined directly by the flow speed in the contraction, so 

that the problem is reduced to finding Cga (r) and X (r), the virtual 
characteristics of the wake. It can be shown that these become 
(Appendix F) 


C»a(r)= a(,j) 


(67) 


X(r) = 



( 68 ) 
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WAVE GENERATOR 




Figure 34 The approximation of the contraction characteristic 


In the limit of a uniform flow, Equations 65 and 66 degenerate to the 
purely dissipative forms and conversely for the zero dissipation situa- 
tion the contraction form alone is recovered. If a contraction charac- 
teristic of the type 


(r') = \u,^(r) - U,„ W] (fff) ^ W 


(69) 


is assumed which corresponds to a linear variation in velocity through 
the contraction, 6b and 66 reduce to 
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These expressions were evaluated and the results shown in Figures 35 
and 36. As was expected, the alleviation of the velocity deficit by 
dissipation is accentuated by the contraction. The wake width, however, 
can either increase or decrease through a contraction depending on the 
relative strengths of the two effects. It should be noted that Equa- 
tions 70 and 71 are the results for a two-dimensional self similar wake 
convected through a one dimensional flow contraction. These expressions 
will be assessed using the results from the test program. 

Consideration of the constancy of circulation in Section 7 leads to the 
relationship between transverse velocities in a vortex upstream and 
downstream of a contraction. It was found that 
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This expression will be used to define the effects of contraction on 
steady or quasi -steady vortical flows. However, these expressions do 
not tell the whole story. A vortex convected through a contraction will 
also experience changes in the axial velocity component. Batchelor (53) 
states that "acceleration of the fluid outside of the vortex leads to a 
change of axial velocity of the fluid within the vortex which is in the 
same direction and of greater magnitude than that outside". This effect 
will be studied in greater detail once test data from the test program 
becomes available. 
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Figure 36 The contraction of a viscous wake on the effect on wake 
velocity deficit. 
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9.3 The Effect of a Screen on a Steady Distortion 

The modification of a steady distortion by a screen is described by 
Equations 41 and 43, which were developed to define the effects of the 
screen on turbulence. Figures 37 and 38 show how the steady velocity 
components are attenuated by honeycomb and gauze type screens in terms 
of their characteristic parameters (length to diameter ratio and so- 
lidity, respectively). 

The constraints on these expressions are that the distortions are small 
compared to the mean flow velocity at the screen. In the static test 
situation, this constraint is probably satisfied by the distortion in 
the streamwise velocity component. However, the transverse velocities 
encountered in a ground vortex may, at the screen be substantial. 

Again, the data from the test program will be used to modify the model 
as required to account for this. 

9.4 Conclusions 

In summary, the effect of contraction and screening on steady or quasi- 
steady distortions is modelled by Equations 70, 71 and 72 for contrac- 
tion and Equations 41 and 43 for screening. 
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DISTORTION VELOCITY RATIO 



Figure 37 The distortion attenuation characteristics of honeycomb. 
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10.0 THE ACOUSTIC TRANSMISSION LOSS OF AN INFLOW 
CONTROL SCREEN (PHASE II TASK E) 

10.1 Introduction 

The inflow control screen mounted upstream of the inlet constitutes an 
obstruction to the acoustic radiation from the engine inlet. The nature 
of this obstruction and its effect on the radiation field are deter- 
mined by the geometrical and structural chacteri sties of the screen 
together with the mean throughflow, and the incident sound field. The 
effect of this obstruction is to produce a reflected field on the 
source side of the screen in addition to the incident field, and on the 
opposite side of the screen, a transmitted field. The interaction of 
the sound field and the screen is very complex and the problem has not, 
in general, been solved. Fortunately, if certain assumptions are made 
that are reasonable for some practical situations, the problem becomes 
tractable. 

The first assumption is that the wavelength of the incident field is 
large compared to the screen aperture size. Secondly, the number of 
apertures per unit area is assumed sufficiently large and uniformly 
distributed so that at some small distance € , on either side of the 
screen, the radiation impedance is constant. (Figure 39) 

The consequence of the first requirement is that the pressure is in 
phase over the aperture and thus the particle velocity of the fluid in 
the aperture is approximately normal to the screen plane. This leads to 
the view of a "locally reacting" screen, the characteristic of which is 
that in the vicinity of the screen for any given condition, the pres- 
sure and particle velocity at a point on one side are solely related to 
those guantities at the corresponding point on the opposite side (Fig- 
ure 39). The second assumption provides a plane on either side of the 
screen containing all of these points and thus the acoustic character- 
istics of the screen can be defined in terms of the radiation impedance 
at the € planes; i.e., two complex numbers. In addition, the locally 
reacting screen assumption, leads to the deduction that the form of the 
transmitted field is the same as that of the incident sound field, its 
amplitude and phase being different. A further reasonable assumption is 
commonly made for transmission through thin screens, namely, that the 
wavelength of the incident field is large compared to the screen thick- 
ness. The screen is thus assumed to be acoustically thin. In this 
situation the pressure is approximately in phase through the screen and 
the fluid moves as a solid body in the aperture. When this occurs, the 
continuity equation provides a relationship between the particle velo- 
cities at the € planes (they are equal in amplitude and phase) and the 
acoustic character of the screen is conventionally described in terms 
of a single quantity, the screen impedance. 
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In view of the wide-spread use of perforated plate as a lining material 
much work has been done on the impedance on this type of structure. The 
dependence of the impedance on the various pertinent parameters, i.e., 
flow, incident sound level and screen geometry has been examined. In- 
variably a single orifice has been examined (54), (55) from which using 
the second and third assumptions the screen impedance is readily de- 
duced. The most commonly considered flow condition has been a grazing 
flow since this is encountered much in practice (56), (57). However, 
Ingard and I sing (54) among others, have made measurements of the im- 
pedance of orifices in the presence of normally incident flow. Their 
measurements of the resistive component of the impedance indicate, 
above a certain flow speed, a linear dependance on that quantity. This 
lends weight to the view that in this region the acoustic pressure loss 
is due to the same mechanism (jetting) as the mean pressure loss and 
can thus be determined from the pressure drop coefficient. At speeds 
approaching and including zero, the viscous dissipation in the acoustic 
boundary layer controls the resistance term. The reactive component, 
they found, is approximately the same as the no- flow value at low 
speeds. It should be noted that the reactive part of the impedance is 
determined by the inertia of the air in the holes, and since the air 
adjacent to the screen at the apertures is also carried along, there is 
an additional contribution to the inertia that must be accounted for. 
These so called end corrections are thus apparently unchanged by a low 
incident throughflow, in contrast to the situation of a grazing flow in 
which the attached air mass is blown away at very low flow speeds (56). 

For honeycomb, the large wavelength to cell depth requirement may or 
may not be satisfied. In the former situation the acoustically thin 
screen approach is directly applicable provided the current screen im- 
pedance is determined. If the wavelength is not large compared to the 
cell depth the flow in the cell may no longer be considered incompress- 
ible. However, the "locally reacting" phenomenon leads to the notion 
that the field propagates through the cells as plane waves and this 
observation allows the problem to be solved. 

On the test stand for engines in the size range of the JT9D and JT15D 
it is anticipated that the screen design will be such that the assump- 
tion of a locally reacting and acoustically homogeneous screen will be 
appropriate. 

10.2 Theory 

For the full scale test situation it is anticipated that the screen 
will be in the acoustic far field of the engine for a substantial part 
of the frequency range of interest. For the JT9D engine this criterion 
is satisfied for frequencies above roughly 1000 Hz if the screening is 
1.5 diameters away. For' the JT15D with the screening located at the 
same number of diameters away, the far field criterion is satisfied 
only above around 3500 Hz. Furthermore in most ICS designs, the screen 
elements are far enough from the source that the incident sound field 
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may be represented by a summation of plane waves and thus, in view of 
the linearity of the acoustic fields being examined, it is sufficient 
to solve the problem for one incident plane wave only. This is done for 
the acoustically thin screen (perforated plate) and also for honeycomb. 
The assumption of locally reacting homogeneous screens is, of course, 
also made. 

10.2.1 The Acoustically Thin Screen (Perforated Plate) 

For a normally incident flow (Figure 40) the convected wave equation is 
satisfied on both the upstream and downstream side of the screen by 
plane waves of form 







(73) 


The corresponding particle velocity field is 
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With, in particular, the magnitude of the component in the x-direction 
being given as. 
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The sign depending on the direction of wave propagation. The acous- 
tically thin screen assumption, which implies that the fluid between 
the € planes is incompressible, allows the normal particle velocities 
on these surfaces to be equated. 
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Figure 39 The assumptioxis for the screen transmission loss analysis. 
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The locally reacting screen produces a relationship between transmitted 
and incident fields so that 



(77) 


The matching of the trace wavelength on the incident € plane for the 
reflected and incident waves together with the constant value of fre- 
quency in these two fields leads to 
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Finally the acoustic pressure balance between the € planes yields 
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Where ^ is the screen impedance 

Solving these equations simultaneously allows expressions for the 
transmission and reflection coefficient to be found. 
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where 
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For the limiting case when M « 1 which corresponds to the situation 
at inflow control structures, >( li. The normally incident wave is 

attenuated most and this occurs when 6^=^0 . Incorporating these two 
conditions into Equation 81 and taking the absolute value of the trans- 
mission coefficient 




. X 


(83) 


where has been replaced by 

Now if the resistance is due to jetting then 
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and from Ingard and Ising the reactive component may be written 


^ ^ (85) 


Those authors note that the end correction S varies from .85d to zero 
as the open area ratio o- goes from high to low solidity. This expres- 
sion becomes increasingly in error as the limits of cr are reached. 
However, in the mid solidity range 6 may be approximated by (l - cr)cLy 
thus Equation 85 becomes 
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substituting the expressions for^ and X. in Equation 83 produces 
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This transmission loss expression is plotted in terms of dB for a 50% 
open area screen in Figure 41. In Figure 42 the loss as a function of 
frequency is shown for the perforated plate used in the JNRP developed 
inflow control screens. The perforated plate having 51% open area, .3cm 
thick and with .47cm diameter holes has significant transmission loss 
at the higher frequencies while the perforated plate having 46% open 
area, ,08cm thick and with .16cm diameter holes has a loss less then 1 
dB at all frequencies of interest. Even though these curves show that 
normal incidence transmission characteristics, the range of incidence 
angles encountered in practice will not materially alter these curves. 

If the criterion for a perforated plate is that the transmission loss 
must be less than 1 dB at all frequencies, and if the maximum frequency 
of interest is 10 kHz, Equation 88 leads to 
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This expression is subject to the constraints of the analysis and the 
range over which the reactance expression is valid. The honeycomb pre- 
sents a different problem for invariably it cannot be considered acous- 
tically thin over the frequency range of interest. 

10.2.2 The Acoustically Thick and Low Solidity Screen (Honeycomb) 

The honeycomb screen problem* while it can be formulated in terms of 
the radiation impedances at its surfaces, is not particularly illumina- 
ted by this approach. The assumption of a locally reacting acoustically 
homogeneous screen results in the sound field propagating through the 
honeycomb cells as plane waves. In addition the low mean flows encoun- 
tered in practice by these screens produce an angle of reflection that 
is equal, to a good approximation, to the angle of incidence (Figure 
43). Viscous effects will be ignored in view of the low pressure drop 
coefficient and a single reflection only will be assumed to occur in 
the honeycomb cells. 

The plane waves in the incident and reflected fields are of the form 
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TRANSMISSION LOSS -dB 



a = 0.51 
t - 0.305 cm 
d = 0.475 cm 


o = 0.46 
t = 0.079 cm 
d = 0.157 cm 


Figure 42 The transmission loss of JNRP developed perforated plate 



with axial components 



(91) 


The transmitted wave while of the same form is displaced by some phase 
that at this point is unknown and hence that pressure field can be rep- 
resented by 







(92) 


where is complex. 

In the honeycomb cells, the waves are of the form 
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Now at the ends of the honeycomb (which correspond very closely to 
the € planes) both acoustic pressure and normal velocity fields must 
match and so at x, = o 


f*z = p/ i-pi 

pr cas ^ - pjf cos 0 = p, 
At X != A 

p,e t p^e 


PrcosS 


(95) 

(96) 


(97) 

(98) 


101 



These four equations are in five unknowns and thus can be solved simul- 
taneously for the transmission and reflection coefficient 

^ cos ^B t sin ^B sin ( 99) 
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In the angle of incidence range encountered by the inflow control 
structures this transmission coefficient is effectively unity for all 
frequencies. As for the perforated plate this observation is qualified 
by the assumptions on which the analysis rests. 

10.3 Conclusions 
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The experience for the acoustic transmission loss of perforated plate 
type and honeycomb type screens that are presented here are 
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for perforated plate 


Ji 0 


TL = !0 

for honeycomb. 


/ - 4 . B sin^ykiJi j 


At the low mean flow speeds to which an inflow control structure is 
exposed the honeycomb is acoustically transparent and perforated plate 
has a transmission loss that is controlled by the reactance of the 
screen. In view of the frequency dependance of the perforated plate 
reactance, significant transmission loss can occur at the higher fre- 
quencies . 
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n.O THE DIRECTIVITY EFFECTS OF SCREEN COMPONENTS 
(PHASE II TASK) 

11.1 Introduction 

Whilo the transmission loss associated with both honeycomb and gauze 
type screens can be confined to 1 dB, the structure supporting the pan- 
els can provide some modification of the radiation field. In fact, any 
discontinuity in the inflow control screen will cause some change in 
the directivity pattern of the source. Ideally, the screens should be 
constructed such that discontinuities do not exist, but within the con- 
straints imposed when operating on a test stand, this is not possible 
although the inflow control structures designed and built under the 
JNRP come pretty close to the ideal in this respect. Other smaller 
scale screens (27, 28, 29, and 31) also exhibit a minimal amount of 
discontinuity in their structure. 

In structures where there are significant discontinuities in the pro- 
file, cau?ed by either support structure or the intersection of honey- 
comb sections, some modification to the radiation field is anticipated. 
The magnitude and extent of the directivity change is dependent on the 
source, so that different sources radiating through the same screen 
will experience different directivity changes. In general, then, an 
engine operating at different conditions will produce radiation fields 
that are modified differently by the screen at each condition. The 
model developed in this section is of a simple source radiating through 
a corner between two adjacent honeycomb panels. The accuracy of this 
model is determined by comparison with measurements made of the phe- 
nomenon. 

11.2 The Model 

The model described here will be used to predict the effects on acous- 
tic directivity patterns of ICS honeycomb panel intersections. This 
model may also be used to determine the effect of structural members on 
the sound field, however, the present discussion is confined to a 
honeycomb corner. 

The source of sound at the test stand is, of course, the engine, and 
since the size of this source is not small compared to the distance to 
the screen, the source cannot be considered to occur at a point. If 
however, this source can be represented by a distribution of simple 
sources and if further the wave equation governs the propagation, then 
the resultant radiation field may be constructed by superposition of 
the radiation fields due to each source. Now, by the time a wave emit- 
ted from one of these sources reaches a corner, its radius of curvature 
is large compared to a typical corner dimension and therefore may be 
considered locally plane. Furthermore, as noted in Section 11, if the 
wavelength of the wave is large compared to the honeycomb cell size, 
the wave propagates through each cell as a plane wave, then reforms on 
the other side as a wave whose direction is unchanged. Figure 44. 
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Consider now a field point that receives significant contributions to 
the fluctuating pressure level through each of the adjacent panels. 
Usually, the source will not lie on the line of symmetry of the corner, 
and so rays passing through the adjacent panels will travel along dif- 
ferent paths, with different path lengths to reach the field point. The 
two rays will consequently interfere with each other after they have 
passed through the honeycomb structure (Figure 44) (N.B. for simplici- 
ty, Figure 44 shows the source and field point lying in a plane normal 
to the corner panels). In general, this is not the case. The nature of 
the interference depends on the source signal for a given geometry and 
for linear radiation, it is sufficient to solve the problem of a point 
monopole radiating at a single frequency. More complex sources can be 
determined by superposition. 


The phenomenon that will be modelled is then the interference effects 
experienced in the far field of a monochromatic source due to the pres- 
ence of a honeycomb corner. The model chosen is shown in Figure 45 
(N.B. again, as in Figure 44, for simplicity, both source and field 
point are shown lying in a plane normal to the corner panels, the model 
is not limited to this case). It is three dimensional. It consists of 
two parallel semi-infinite and totally reflecting planes, on the oppos- 
ing sides of which are two simple sources of the same frequency. The 
distance between the planes is related to the blockage of the corner. 
The phase relationships of the two sources is determined by the path 
length difference of the waves passing through the adjacent panels. The 
distance of the sources from the plane boundary is the distance from 
the source to the corner. 

The model then, exhibits all the important characteristics of the phe- 
nomenon. 


To find the solution to this problem, a look at the diffraction pattern 
for a single source and plate is beneficial (Figure 46) because the two 
plate solution may be constructed from the single plate solution, by 
super-position, in view of the linearity of the problem. 

The solution is presented by Bowman, Senior .and Uslenghi in (58) and is 
given for the velocity potential at the field point (Jr,$). 


* (R, *) = 


„ Hf (/i"- i- '-kR) J 


(101) 


where 

m = ± 


"h VS. for 0<^^TT 
~ VS for 


This expression may be simplified (Appendix F) by requiring that 


2 « /? , JlR» J 
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(The validity of these and subsequent assumptions in this section will 
be examined in Section 11.4). 


A simple expression for the velocity potential can be obtained in terms 
of Fresnel integrals and is 


!■< -^)| |(x 


( 102 ) 


The Fresnel integrals C and S , are functions of m , which for a 
fixed value of R. and ^ means that they are functions only of the 
field point angular position 0 and the coordinate along the 
corner Z . The expressions for them are 


( 7 * / ^ £U>S T ^ ) c/ T 


VF 


(103) 




Now noting that the intensity in the far field is proportional to the 
square of the amplitude of the velocity potential (i.e.0'0'* ), the 
intensity can be written as shown in Equation 104, and may be plotted 
as a function of 0 for a given R, and ? . A typical shape is 
shown in Figure 47. ‘In general the shadow is not sharp edged and the 
boundary is characterized by a diffraction pattern with interference 
fringes on the source side of the plane. This shape is due to the 
Fresnel integrals. 


Intensity 





(104) 


It is noted, that this form for the intensity is identical to that pro- 
duced for the case of Fresnel diffraction around a half plane, (59). 
Apparently the constitutes a definition of Fresnel 

diffraction for this geometry. 


Using the solution to the single plane problem the solution for the 
model. Figure 48, can be found, by superposition, and in terms of the 
velocity potential is 






jCi- 


(105) 


-t 
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The subscripts refer to the upper and lower sources, and t is 
the phase difference between the sources. 


There are three distinct regions of the field - firstly, where /r>uis -f-i/a 
and mji is -ve , corresponding to field points that are exposed to the 
incident fields of the upper source but in the shadow of the lower. The 
second regime is the converse, while the third consists of those posi- 
tions in the geometric shadow of both (Figure 48). The intensity can 
thus be written, (within a constant) as 



(106) 




i~ Ou) ^ Cj) Hii-S^Xi +St) 




f 






The first two terms are the separate contributions from each source, 
while the third and fourth are due to the interference between them. 

In This expression all distances appear as a product with the wave num- 
oer K . For ease of computation the field point is referenced to an 
origin halfway betwen the two sources and the relevant relationships 
are shown in Figure 48. 

At this point, it is necessary to relate the parameters of the model to 
the geometry of the honeycomb corner. The plate separation, h and 
source phase difference f , may be determined from the geometry of Fig- 
ure 4b in terms of the corner angle, /3 , honeycomb depth, X , and 
incidence angle 0 so that 

h ~ 2. cos cos ^ 

“ Zl. cos ^ (107) 

There is an assumption implicit in the above expression for the phase 
factor; it is that ray tracing is valid for determining this factor. To 
Quantify the range of validity, it is noted that geometric acoustic 
techniques are valid when ■fia»l where a is a characteristic dimension 
normal to the propagation direction. This is also the region in which 
scattering is referred to as diffraction. Since as noted previously, 
the model is a diffraction model, when the condition that^A be large 
IS not satisfied, neither the phase calculation nor the use of the 
model is valid. 


The intensity in Equation 106 is divided by a reference intensity and 
converted to dB. This process is described in the following section. 
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11.3 Comparison Between Experiment and Theory 

As part of Phase II, Task F, measurements were made to determine the 
directivity effect of a honeycomb corner with a view to assessing the 
theoretical model. This test program was performed in an anechoic cham- 
ber at X-206 stand. Figure 49. 

Briefly, a honeycomb corner was set up in the anechoic chamber and ex- 
posed to a driver .91cm feet away (Figure 50). An array of microphones 
at a radius of 15 feet was used fo measure the transmitted sound field. 
The incident angle, Q , could be varied between 0°, 5°, and 10°. 

The honeycomb was 7.62cm inches thick with an included angle of 135° 
between panels. Cell size was .95cm in diameter. The source was driven 
at nominally, 2kHz, 5kHz, and 9kHz and was not omnidirectional in the 
angle range of interest. The corner was covered by tape in order to 
prevent sound leaking through it. 

Two types of tape were used; a green cloth tape and an aluminum tape. 

No substantial difference was observed between data gathered from the 
two situations, except for the 2kHz case at 5° incidence. Here it is 
conjectured that an inadvertant configuration change was made in the 
aluminum tape case (e.g., moving of absorbing material in the anechoic 
chamber). 

The measured directivity patterns are shown in Figures 51 through 59. 
These were obtained by subtracting the levels obtained with the corner 
in place from those obtained without it. In general, the directivity 
changes increase with the frequency, indicating that the corner is more 
disruptive to the radiation field as the incident sound wavelengths 
decrease. In addition, there is, at the two higher frequencies, a local 
maximum in sound pressure level directly behind the corner as viewed 
from the source, at zero incidence there is no shadow zone behind the 
corner, the directivity patterns being seemingly arbitrary to the in- 
tuitive eye. These patterns are compared with those obtained from the 
model in these same figures. 

To obtain directivity patterns from the model it was necessary to de- 
termine a reference level. This was done by observing that the inter- 
ference occurred predominantly at a range of angles in the vicinity of 
the source - corner axis, and outside of this range the intensity lev- 
els were substantially constant. It was this constant level that was 
used as the reference (or corner removed) level . 

In the table of Figure 60, two quantities are shown that provide an 
assessment of the model. The first is the mean difference between the 
measured and predicted directivity patterns. The second is the standard 
deviation of the measured data from the theoretical directivity pat- 
tern. These quantities are computed for each configuration and are used 
together with the comparisons of measured and theoretical predictions, 
shown in Figures 51 through 59 to provide the following assessment of 
the model . 
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P&WA X-206 stand anechoic chamber 



Figure 49 A schematic diagram of the anechoic chamber. 
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Figure 51 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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Figure 52 Diffraction at honeycomb corner comparison of measurement 
and prediction. 


115 


AdB REDUCTION 


f = 2000 (1995) Hz ^^ = 3' 


O fi A 


8 6 ® 


O EXPERIMENTAL 
(GREEN TAPE) 

A EXPERIMENTAL 
(ALUM TAPE) 

-- THEORETICAL 


O /!- 
/ 

/ 

AN( 


T110° 


ANGLE 


ANGLE 


Figure 53 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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Figure 54 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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Figure 55 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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Figure 56 Diffraction at tx>neycomb corner comparison of measurement 
and prediction. 
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gure 57 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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Figure 58 Diffraction at honeycomb corner comparison of measurement 
and prediction. 
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ANGLE OF 
INCIDENCE 
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0.45 
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0.79 

1,138 
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FIGURE 11.3.12 VALUES OF (i) THE MEAN DIFFERENCE BETWEEN MEASUREMENT AND 

THEORETICAL PREDICTION 

(ii) THE STANDARD DEVIATION OF THE MEASUREMENT FROM 
THE THEORETICAL PREDICTION IN dB. 


Figure 60 The values of (i) the mean difference between measurement 
and prediction, (ii) the standard deviation of the 
measurement from the prediction in dB. 
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At 2kHz, all cases indicate that the levels are predicted well in an 
absolute sense and that the shape is well reproduced. The 5° inci- 
dence case is not predicted as well on either count as the other two. 
This is due to 3 data points around 130° taken for the aluminum tape 
configuration and these points are in conflict with the green tape mea- 
surements, which indicates that there were other differences between 
the configurations. 


At 5kHz, both absolute levels and shape are very well predicted. In 
this case, it is apparent that the model is very impressive. 

At 9kHz, a superficial glance indicates that the shapes are reasonably 
well predicted with absolute levels not very impressive, however, there 
are^some alleviating considerations. Firstly, at this frequency small 
variations in angular position result in large changes in level, 
pe position accuracy of the microphone has been estimated at 1/2°. 

Such a variation would cause large variations in the two parameters of 
Figure 60 particularly for the 10° incidence case. Thus at this high 
frequency case, errors can be attributed to measurement position error 
as well as model error, and discrimination between them is difficult. 

Neglecting all alleviating considerations, however, upon averaging over 
all tests, 74% of the data was within 1 dB of the predictions. The IdB 
figure is chosen since this is the error bound of measured SPL's. In 
addition, the mean difference obtained by averaging over all tests is 
about .38 dB and the corresponding standard deviation from the model 
prediction is about .9 dB. 

11.4 Validity of the Model Assumptions 

(i) 2^o»R- 

For the full scale inflow control structure, fitted on a 
JT9D, at a typical stand, (the source to screen corner 
distance) is about 4.57m. (the microphone array radius) is 
4.57m. The assumption is thus marginally valid. 

For the experimental situation ^owas .91m, and/^was 4.57m. 
The assumption is not valid in this case. 

(ii) Acl»I 

For the full scale measurement situation, at a frequency of 
2kHz (approximately B.P.F. at approach for JT9D), 1700. 

The assumption is definitely valid. 

For the experimental situation, at 2kHz, 5 : 170. Again the 
assumption is valid. 
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(iii) -^^L»/ 

For the full scale inflow control structure, cl is about 
7.62cm. and at 2kHz, Acl'^3. This was also the lowest 
value for the experimental situation. This assumption is 
not valid at this frequency, but at frequencies above 5 kHz 
it is probably valid. 

11.5 Concluding Statements 

In spite of the violation of two model assumptions by the experimental 
arrangement, the theory produces predictions that are quite close to 
the measured directivity patterns. In addition, the full test situation 
will tend to meet one of these assumptions better. Consequently, this 
model can be used with some confidence in the full scale case. 

The application of this model to a distributed source can be implemen- 
ted by summing the contributions to the far field from each individual 
source region. A determination of the radiation field due to a line of 
sources is currently being made. 
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12.0 IMPLICATIONS OF ANALYTICAL STUDIES TO ICS DESIGN 

Based on the analytical studies reported herein, the following general 
observations can be made to guide future ICS designs. 

1. It is becoming increasingly apparent that the inflow distor- 
tion field is not solely due to direct contraction of atmos- 
pheric turbulence but also receives contributions from flow 
over stand structure, ground effects, etc. The consequence of 
this is that, to produce an optimum screen design for a given 
engine, the immediate structural environment must be accoun- 
ted for. This obviously is going to vary from stand to stand. 
The most reasonable design course, given the variation in 
stands seems to be to specify a worst case and design for 
that. It should be noted that honeycomb and perforated plate 
suppress steady as well as unsteady velocity distortions. 

2. In existing inflow control structures various combinations of 
honeycomb and perforated plate are used. The rationale for 
using honeycomb is quite easy to understand since transverse 
velocity components are amplified substantially by the sub- 
sequent contraction. The use of perforated plate which sup- 
presses preferenti ally the streamwise components is somewhat 
harder to justify since the contraction strongly attenuates 
this component. However, the rapid initial rate of return to 
isotrophy downstream of the screen will cause the transverse 
velocity component to increase in value and thus be subse- 
quently amplifiable by the contraction. The suppression of 
streamwise distortions appears then desirable in order to 
prevent this transfer of energy back into the transverse dis- 
tortion component. 

3. In order to render honeycomb transmission losses negligible, 
it is necessary to design for small through-flow velocities, 
on the order of 4.63 - 9.25m/s. Designs of this type require 
that the honeycomb structure be located outside the inlet, 
and in a region where inflow velocities are low. 

4. Perforated plate transmission losses can be minimized by in- 
voking Item 3 and requiring that design criterion of Equa- 
tions 89 be met. 

5. Directivity effects can be minimized by designing ICS's with 
a minimum of structure and discontinuities in the honeycomb 
and perforated plate surfaces. 
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13.0 CONCLUSIONS 


1. A transfer function was developed that relates blade mounted 
transducer pressure measurements to inflow normal velocity. 
This function is defined in Figure 13 and will be applied to 
define inflight and static fan velocity distortion field. 

2. Observations of PWA/Boeing Joint Noise Program Blade Mounted 
Transducer Data show that during static engine testing sig- 
nificant inflow distortions can result not only from atmos- 
pheric turbulence, but from the ground stand structure as 
well. Furthermore, the ICS's designed by PWA and Boeing for 
application to the JT9D engine were reasonably effective in 
removing these distortion types. These findings are discussed 
in Section 6 of this report. 

3. Existing theoretical studies on the effects of contraction 
and screening on flow distortion have permitted formulation 
of analytical models that account for the effects of contrac- 
tion and screening on steady and unsteady flow distortions. 

The models for the effects of contraction on unsteady and 
steady distortions are given in Equations 20, 21, 22, and 70, 
71, 72, respectively. The models for the effects of screening 
on unsteady and steady distortions are given by Equations 45, 
46, 48, 49, and 41, 43, respectively. 

4. Models accounting for the effects of honeycomb and screen on 
acoustic transmission loss were developed and are given by 
Equations 100 and 88. 

5. A model accounting for the effects of honeycomb and structure 
on acoustic directivity patterns was developed and shown to 
agree well with experimental data. This model is given by 
Equation 106. 

6. Based on the studies carried out in Phases I and II of this 
contract, guidelines for future ICS designs were identified. 
These guidelines are discussed in Section 12 of this report. 
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APPENDIX A 

The Power Spectral Density of the Unsteady Upwash Velocity Component of 
the Inflow Field 

To determine the blade mounted transducer (BMT) response function it is 
necessary in part to relate the upwash velocity power spectral density 
to the measurements made by a hot film probe. This probe measures the 
streamwise and the circumferential component of the unsteady velocity 
field. The relationships between these various velocity components is 
shown in Figure 4. 

By considering the rotation of the coordinate system it is seen that 


t = - u, (t-) cos -fUi (t) sin y 


( 108 ) 


Taking the Fourier transform of 108 yields 


(w) (j) Cos y vS5 Cw) sin Y 


(109) 


The power spectral density is proportional to J^(co) 3^ (cS) so that the 
relationships between the power spectral densities of the various ve- 
locity components is 


(u) - F] (<j) cos ^ y i" Fz Cw) s/n ^ If - 


( 110 ) 


Kobson (60) has considered the general problem of the effect of angular 
rotation on spectral quantities and derives an expression corresponding 
to no. There is, however, a misprint in that expression. 

Estimates of the various terms in Equation 110 on an order of magnitude 
basis for the case of static operation of a JT9D at an approach power 
setting allows some simplification of the expression. 

According to hot film measurement and analysis, in the tip region not 
in the boundary layer, 

F^M-o \iof,m\ ^0 \.oi\ 

' • ' * ( 111 ) 
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the magnitude of the real part of the cross spectral density func- 
tion must be less than the total magnitude of the function or. 


Rj. (-c^W||< \f,2U 


012 ) 


but from Equation 111 




(113) 


Therefore, 



(114) 


Now in addition, in the vicinity of the tip the blade relative flow 
angle, ^ ^ 3Q0. 

Thus the order of magnitude of the terms in Equation 110 become, for 
the static JT9D case under consideration; 

F^(J) = ,7s ol^/o io 


F„(J) = O ^3.25 FiM]^ 10^.3F,M^ 

Therefore the term containing the real part of the cross spectral den- 
sity is small compared to the sum of the other two terms and thus may 
be neglected. For this case, then Equation 110 may be simplified to. 


(w) — (J) cos ^ H i" sin ^ if' 


(117) 
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APPENDIX B 


Criteria for the Linearization of Equation 16 and the Applicability of 
that Equation to the JT9D Inflight and Static Operation 

As stated previously the equation of motion for an incompressible, con- 
stant viscosity unsteady flow with zero mean vorticity is from Equation 
16. ^ 


Dt 




(118) 


Taylor linearized this equation by neglecting the last two terms of 
Equation 118 thereby recovering Cauchy's equations. The justification 
for neglecting these terms is that of considering only a rapid distor- 
tion. Batchelor argues that if the contraction is sufficiently rapid, 
the straining of the fluid by the mean flow is the dominant effect. He 
formulates the criterion for the rapid distortion as 


i=r- 



(119) 


The change in the unsteady kinetic energy of the flow is that due to 
the turbulence inertial effect and the viscous decay Ui average is 
the mean particle velocity through the contraction. Equation 119 
states that the rate of change of the turbulent kinetic energy due to 
the two above mentioned effects is muclijess than the contraction rate. 
Batchelor proceeds to approximate with an expression empiri- 

cally derived for turbulence in the initial period of delay and in the 
absence of mean flow gradients. 



( 120 ) 


Where -t is a turbulent scale and . 

Substituting this expression in Equation 119 and assigning characteris- 

kinetic energy and turbulence length scale produce 
the Batchelor's criterion for neglecting the last two terms in Equation 


In terms of the initial mean flow velocity and contraction 


IAr!±\ 

D \ 2 ) 


(122) 


It is noted that this inequality is always satisfied for sufficiently 
large contraction rate, . Furthermore, on consideration of this 
inequality at any point x , from the beginning of the contraction 


-e. JMt! 

U,„ X, z 


(123) 


It is immediately apparent that this inequality is satisfied for suf- 
ficiently small -X/ . The criterion for linearization in 122 depends on 
the empirical relationship 120 whose appl icabil ity to turbulence in the 
presence of mean flow gradients must be questioned. It is suggested as 
an alternative that the approach used by Ribner and Tucker (37) may be 
more satisfactory. Those authors generate the criteria directly from 
Equation 118. Thus, for neglect of the turbulence inertial and viscous 
decay terms in Equation 118 it is necessary that 


and 






(124) 
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(125) 


Considering first the turbulence inertial criteria 124, Ribner and 
Tucker suggest ^hat the dominant mean flow gradients in a typical con- 
traction are and and these are of comparable 

magnitude. If it is assumed that the (fther gradients are zero, in- 
equalities 124 reduce to 
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(126) 
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To determine more useful forms of the inequalities 126 the gradients 
will be estimated using average values. This 126 becomes 



(127) 


’ "^ 2 » '('3 are Cartesian length scales. In an axisymmetric 
flow the last two inequalities become the same. 


Now the results from the simple contraction model are reasonably close 
to the mean square velocity ratios determined by the Batchelor and 
Proudman, Kibner and Tucker solutions to the problem. Consequently 
these simple model results may be considered approximations to the lin- 
earized theory and as such when substituted in 127 will provide cri- 
teria for the neglect of the turbulence inertial terms. Before this 
substitution is made, however, it is noted that the inequalities in 127 
must be satisfied everywhere in the contraction in order for the lin- 
earization of Taylor to be valid. Therefore, the inequalities will be 
established at the beginning and end of the contraction and the most 
constraining criterion will be the minimum requirement for neglect of 
the turbulence inertial terms. If then the turbulence field is assumed 
initially isotropic, at the beginning of the contraction 

_ .T* — 77 ^ 


and if isotropic turbulence is characterized by spherical eddies then 
and thus inequalities in 127 reduce the 

Ml . M. a D 

^ Ax, ■ T> ^ 

(128) 
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where J. i = contraction ratio 

and B = contraction length 

The A subscript denotes the pre-contraction station. At the end of the 
contraction, using the simple model results in 

f"T.) = 


<^f» 

-U., uf/i 



11 

11 



“ ^ •^SA - J , 






= = 
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Substitution in 127 yields 






U/-/) 

D 31, 
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Comparing the criteria(128, 129a and 129b) it is noted that if 129b is 
satisfied, all are satisfied and so this inequal ity must be satisfied 
if the turbulence inertial terms are to be neglected through the con- 
traction. 

Kibner and Tucker developed the criterion for neglect of the viscous 
decay term in Equation 118 by simply requiring that 


<i>X 




(130) 


and then considering a single shear wave they obtained the result that 


v'A. 


Ml 

^X/ 


(131) 


where is the absolute value of the wave number vector. 


If the mean flow gradient in Equation 131 is approximated, that in- 
equality becomes 



(132) 


It IS the opinion of Kibner and Tucker that the major part of the spec- 
trum will satisfy this condition. It should be noted that the viscous 
neglect criterion 132 is conditional on the turbulence inertial term in 
Equation 118 being negligible. 
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Application of the turbulence inertial effect neglect criterion 129b to 
a JT9D operating statically and in flight provides information on the 
applicability of the linearized theories to these situation. Both sides 
of the inequality are shown in the table below for both operating en- 
vironments. The values of initial turbulence intensity are determined 
from the atmospheric model presented in the Phase I report of this con- 
tract. 



Contraction 

Contraction 

Initial 

Initial 



Ratio 

Length 

Scale 




D 


D 3X, 







Inf 1 ight 

1.5 

3m 

.43m 

.0016 

.016 

Statically 

50 

15m 

3.54m 

.2 

.077 


The inequalities are satisfied for the flight case but not for static 
operation. Therefore the linear solution approximated by the eddy model 
is valid for the in-flight contraction phenomenon. In static operation 
however, turbulence inertial effects are not negligible and consequent- 
ly the tendency to isotropy will alleviate the distortion of the turbu- 
lence due to the flow contraction. In this latter case however, it is 
possible that a reasonable estimate of the turbulent kinetic energy 
change through the contraction is provided by the eddy model if a decay 
correction is incorporated (43). 

The test program of Task F will supply information on the effect of 
non-linear contractions on turbulence, since the linearizing require- 
ment is not satisfied for the small test facility. 
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APPENDIX C 


The Effect of a Contraction on Turbulence Integral Length Scales Using 
the Linear Theory 

Using the results of Batchelor and Proudman (38) and their notation 


0// (^* 0t/ 


(133) 


where 

is the three dimensional spectrum tensor for the ith velocity 
component. 

The prime indicates the pre contraction location 

Ai>Ki the pre and post contraction wavenumber vectors. 

are the diagonal components of the strain tensor. 

with corresponding expressions for ^22 , and (Hi~ O), 

Integrating over and Kjf produces the relationships between the 
corresponding upstream and downstream one dimensional spectra 

K,= 0 . 


p;u,^o) 


and noting that the scale of the ith velocity component in the ^direc- 
tion. • 





(135) 


then 


hi 1 jiL 

L'n ■ < 


(136) 
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and for the other components 

^ it ^2 ^i. ^33 ^3 


(137) 


for an axisymmetric contraction such that then > And 

then 


_ / • ^X£ _ j ^33 3 y-^t 

Ln ^tt ^2 


(138) 


where 



For an axisymmetric contraction it is also true that 


0g2 (^j f ~ (139) 


and consequently 


,£l i J=£L ^ djL. . (140) 

^2/ <^Z ^3! ^3 


These relationships are all independent of the initial turbulence 
field, however, for any initially isotropic field Vf * and 
Batchelors' expressions for and v*|x^s|x^ may be used in 

Equations 138 and 140. For the latter case the ratios are shown in 
Figure 19. 
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APPENDIX D 


The Effect of Screens on Turbulence Integral Length Scales Using the 
Linear Theory ^ 

As in the contraction effect, certain length scales of an isotropic 
turbulence field are changed by the screen in a manner that does not 
depend on the particular form of the field. 

Batchelor's formulation of the solution(43)a little manipulation 
shows that 


% 

(0. 

K,, 

K3) 

* (Hi^o) $// (Pj /{tf 

(141) 

03 * 
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Ks, 


* “ ip, Hz >^s) 

(142) 
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- ^'sziO.Hz.K,) 

(143) 

032 

C/C, 

,0, 

K,) = 

■■ Oi,p,Kj) 

(144) 

0,, 

c/c„ 


.0) = 


(145) 

where cTV^ * 

(Kr 

^0) - 

/ /-a — 0 K 

l-l-afK 



0^'^* is the velocity spectrum tensor and 

the prime denotes the downstream location 
i s the wave number vector 
a is the flow angle ratio 

The coefficients in each case being independent of Kx » and thus inte- 
gration over the non zero wavenumber components produces the relation- 
ships between the one dimensional spectra value at zero. 


f’, % ' o') = - d) p, '% = o) 


(146) 
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Pz (K, =o) = 

P2(K,^0) 

(147) 

Pf iKi-0) = 

P,'%-0) 

(148) 

Pz (Kz’O) - 

PziKz-’O) 

(149) 

Pj%--0) = 

II 

(160) 


and again noting that the scale of the ith velocity component in the ^ 
direction ^ 



Px(J<z-0) 

1 (151) 
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For an initially isotropic field ju^sfx^and Taylor and Batchelor's ex- 
pressions Equations 32 and 33 for and may be used in 

Equations 152 through 156 above. It should be noted that the other 
length scale ratios are dependent in form on the initial turbulence 
spectrum and general statements cannot be made about them. 
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APPENDIX E 

A Comparison Between Measured Honeycomb Loss Coefficients and Values 
Determined Using Blasius' Formula for Turbulent Pipe Flow 

The Blasius formula for turbulent flow in smooth pipes is valid for di- 
ameter Reynolds Number less than 10^ and greater than 3.103. It may 
be written in the form of a loss coefficient. 


K = 




1 

d 


(157) 


where /2 a. is the diameter Reynolds Number 
and is the pipe length to diameter ratio 

The only adaptation necessary to apply this expression to a honeycomb 
is to determine the velocity of the flow in the cell Uc . This may be 
written as U/ao incompressible flow, where 5 is the 

solidity of the honeycomb. 

The measurements conducted by Paterson were made for flow speeds, using 
.95cm diameter honeycomb with a length of 7.62cm. The characteristics 
of this honeycomb were -6/d = 8 and s = 0.05. 


The results of the comparison are shown in the table below. 



UfC 

Rjl 


K 

calcul ated 

mea^red 

16 

16.8 

3350 

7.61 

.333 

.365 

30 

31 .6 

6302 

8.91 

.284 

.295 

50 

52.6 

10490 

10.12 

.25 

.26 

68 

71 .6 

14279 

10.93 

.232 

.255 


Table 1 Comparison of Measured and Calculated Honeycomb Loss Factors 

The turbulent pipe flow loss coefficient provides reasonably accurate 
estimates of the zero incidence angel pressure drop coefficient for 
this honeycomb. The pressure drop coefficient for this honeycomb ap- 
pears to be dominated by skin friction, however, for higher solidity 
honeycombs other components of the form drag, e.g., base drag, may be 
important. 
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APPENDIX F 


The Virtual Characteristics of the Dissipative Wake in a Flow Contrac- 
ti on 

The contraction of the wake is approximated by a series of steps in 
which alternately the wake is dissipated in a uniform flow and then 
contracted instantaneously. In each of the dissipating regions the wake 
can be considered to evolve as if it were generated, for each region, 
by a rod of dimeter a., drag coefficient Cg located at a distance X . 
These values will be different for each dissipative region and are thus 
termed here as "virtual characteristics". In each dissipative region 
the product C^a will be invariant and will change instantaneously 
across each contraction. Figure 61. 

The assumption of wake characteristics obeying the Equations 64 allows 
the observation that 


AU,-«- = Cpa 


(168) 


at any point in the flow contraction. Thus across the first instantan- 
eous contraction 


~ A U, (r-g) ^ ^ U,^ (r^ 

C.iO.<r„) Aiv,Cr.) ' ' 

by definition ^/,(the contraction ratio of Step 1) and from 

the approximate contraction theory 

AU,(rJ ^ ^ j_ 

l^Ui(r^) Cj (160) 

CDO-^re) _ .-'/K 

but Cp CL 

Cpair^) ^ CpO. {r^) 

Cjj C, ^ C.pCL 
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and after n contractions in a length (r- r^,) 

.C^ O'*"* i^paCO <'«') 

Now letting n-»oo 

CoQ^ir) * {c, CnY^* a.(n,) 

. -avi 

= i,(r) Cpa(n>) 

where U,ir) is the overall flow contraction ratio to the point r . 

The virtual source distance varies through the dissipative region and 
changes instantaneously at the contraction Figure 61. 

Again from 64 


_ x_ 

^u, ' u,^ 


(163) 


which yields the relationship of the virtual source distance across the 
first contraction 


X(rA) 


= c 


3l2 


(164) 


but 




thus 


(X(r^)-f- A ) 


(165) 


and after n contractions 
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Writing in terms of the mean velocities 



n-/ 


or 


X(r) = Xft-) X(rJ-i- 



( U^(r) 



O- 


(166) 


Let 


Equations 162 and 167 contain the expressions for the virtual wake 
characteristics that are used in Section 9. 
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APPENDIX G 


Simplification of the Velocity Potential Due to a Source on a Half Plane 

Bowman, Senior and Uslenghi(57) present the solution for the non di- 
mensional velocity potential for the hardwall case as 


* (H, *) * 



( 168 ) 


where 

nt — ^ 6 ^/ 


See Figure 46 for the geometry 

For the farfield regime -AR^I , and so 

Now from Morse and Ingard (61), for large arguments the Hankel function 
may be approximated as 


|vj^- ^|j 

(169) 

= Ziyf^ e.Kp 

(1?0) 

Now putting 

W = ZZ y/i e<p {RR- j 


and splitting the integration in two, leads to 


ri 

T, X, 

(171) 

Examining I, ; note that „fn 

= 


Therefore if it is assumed that R then in the range of in- 
tegration, the denominator of the integrand may be considered indepen- 
dent of /4 . 
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Thus 


(172) 


-i -AHy/T dcRVI 


Evaluating this integral in terms of Fresnel integrals, see Equation 
103, yields 


for 

7 - « 

f " 2A^ 


(173) 

and for m-v£. 

7 - _ 

2AR 


(174) 

Then for m -tvt. 

^ ^ 1^/ ^ 

/j-.j 

(175) 

and for m-v& 



(176) 


where \nn\ is used in the Fresnel integrals. Now note that the sum of 
these fields produces the same field as that generated by a simple 
source of twice the strength 


1 .e. 


mi- nt m- y* 


2a 


iAd 


itd 


(177) 


or 




4a 

which gives ^ — 

= m 

Substituting 178 and 175 or 176 produces the expression for 


(178) 


K, s 




riu I/O or I/O proauces une expression to 
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velocity potentiel for a simple source on a half plane 
IS the building block for the honeycomb corner diffraction model of 
Secti on 1 1 . 
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